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ABSTRACT 

We report on the discovery of HAT-P-llb, the smallest radius transiting extrasolar planet (TEP) 
discovered from the ground, and the first hot Neptune discovered to date by transit searches. HAT- 
P-llb orbits the bright (V=9.587) and metal rich ([Fe/H] = +0.31 ±0.05) K4 dwarf star GSC 03561- 
02092 with P = 4.8878162 ±0.0000071 days and produces a transit signal with depth of 4.2mmag; the 
shallowest found by transit searches that is due to a confirmed planet. We present a global analysis of 
the available photometric and radial- velocity data that result in stellar and planetary parameters, with 
simultaneous treatment of systematic variations. The planet, like its near-twin GJ 436b, is somewhat 
larger than Neptune (17M ffi , 3.8i? ffi ) both in mass M p = 0.081 ± 0.009 Mj (25.8 ± 2.9 M e ) and radius 
Rp = 0.422 ± 0.014 R 3 (4.73±0.16i?©). HAT-P-llb orbits in an eccentric orbit with e = 0.198±0.046 
and u> = 355.2 ± 17.3°, causing a reflex motion of its parent star with amplitude 11.6 ± 1.2ms -1 , a 
challenging detection due to the high level of chromospheric activity of the parent star. Our ephemeris 
for the transit events is T c = 2454605.89132 ± 0.00032 (BJD), with duration 0.0957 ± 0.0012d, and 
secondary eclipse epoch of 2454608.96 ± 0.15 d (BJD). The basic stellar parameters of the host star 



are M* = 0.809 



+0.020 
-0.027 



M , i?* = 0.752 ± 0.021 R Q and T cS ± = 4780 ± 50 K. Importantly, HAT- 



P-ll will lie on one of the detectors of the forthcoming Kepler mission; this should make possible 
fruitful investigations of the detailed physical characteristic of both the planet and its parent star at 
unprecedented precision. We discuss an interesting constraint on the eccentricity of the system by the 
transit light curve and stellar parameters. This will be particularly useful for eccentric TEPs with 
low amplitude RV variations in Kepler's field. We also present a blend analysis, that for the first time 
treats the case of a blended transiting hot Jupiter mimicing a transiting hot Neptune, and proves that 
HAT-P-llb is not such a blend. 

Subject headings: planetary systems — stars: individual (HAT-P-11, GSC 03561-02092) techniques: 
spectroscopic, photometric 



1. INTRODUCTION 

Transiting extrasolar planets (TEPs) are uniquely 
valuable for understanding the nature of planetary bod- 
ies, because the transit light curve, combined with pre- 
cise radial velocity (RV) measurements of the reflex mo- 
tion of the parent star, yield unambiguous information 
on the true mass and radius of the planet, assuming that 
the stellar mass and radius are known. By inference, 
it is then possible to investigate the internal structure of 
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these planets, as has been done by several teams trying to 
formulate and match theories to the observed bulk prop- 
erties of known TEPs (e.g. Baraffe et al. 2008; Fortney 
et al. 2007; Burrows et al. 2007; Seager et al. 2007, and 
references therein). The transits across the face of the 
star enable a plethora of scientific follow-up opportuni- 
ties, such as detection of the atmospheres of these plan- 
ets via transmission spectroscopy (Charbonneau et al. 
2002), measurement of the stellar spin axis versus plane- 
tary orbit (Winn et al. 2005; Johnson et al. 2008a) via the 
Rossiter-McLaughlin effect (Rossiter 1924; McLaughlin 
1924), or measurement of their equilibrium temperature 
while they are occulted by their central stars (Charbon- 
neau et al. 2005). 

Photometric searches for TEPs have published some 
50 such objects over the past 8 years 12 , most of these 
with masses and radii in excess of that of Jupiter. Pre- 
viously the smallest mass TEP discovered by the tran- 
sit search method was HAT-P-lb with M = 0.52Mj 
and R — 1.22i?j (Bakos et al. 2007a), recently super- 
seded by the discovery of WASP-ll/HAT-P-10b with 
M = 0.46Afj (West et al. 2008; Bakos et al. 2008). The 
smallest radius planet detected by ground-based tran- 
sit searches was HAT-P-3b (Torres et al. 2007) with 
R = 0.89i?j, and the smallest radius planet from the 

12 http:/ /www.exoplanet.eu/catalog-transit.php 
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space is Corot-7b (Leger et al. 2009). 

In the meantime radial velocity surveys have been 
reaching down to Neptune-mass planets, thanks to high 
precision and high signal-to-noise spectrographs that de- 
liver radial velocities at the ms _1 level over an extended 
time, such as HARPS (Mayor et al. 2003) on the ESO 
3.6 m telescope, or HIRES on Keck (Vogt et al. 1994). It 
was a major advance when Santos et al. (2004) discovered 
the MpSmi = 14 M® planet around /i Ara, McArthur 
et al. (2004) detected a Neptune-mass planet around 
55 Cnc, and Butler et al. (2004) found a <~ 21M® mass 
planet around GJ 436. These were followed by further 
exo-Neptune discoveries, such as the three Neptune plan- 
etary system around HD 69830 found by Lovis et al. 

(2006) . 

Recently the detection threshold of RV searches has 
reached even below that of super-Earths (M < 10M®). 
Rivera et al. (2005) found a ~ 7.5 M® super-Earth or- 
biting the nearby M dwarf GJ 876. Udry et al. (2007) 
found a 5M® and an 8M® mass planet in a triple 
planetary system around GJ 581. Finally, Mayor et 
al. (2008) discovered a triple super-Earth system with 
4.2, 6.9, and 9.2 earth masses around HD 40307. Al- 
together, as of writing, some 20 objects with minimum 
mass M p sin i < 0.1 Mj = 31.8 M® have been detected 
by the RV technique. 

Radial velocity detections are routinely checked 
for transit events by the discovery teams, or 
by the transitsearch.org collaboration of ama- 
teur/professional astronomers (Seagroves et al. 2003). 
A few successful detections have been reported 13 . One 
such case is HD 189733b, a 1.13 Mj planet around a K 
dwarf, discovered by Bouchy et al. (2005) via the RV 
method, and confirmed to transit by the same team 
via the Rossiter-McLaughlin effect, and then via follow- 
up photometric observations. Another example is the 
21.2 day period eccentric planet HD 17156b found by Fis- 
cher et al. (2007), with transits detected by Barbieri et al. 

(2007) through transitsearch.org. A third example, 
HD 149026b, is a transition object between Jupiter-mass 
and Neptune-mass planets, in the sense that it is a hot 
Saturn with M p — 0.36 Mj (or about 1.2 times the mass 
of Saturn), and R p = 0.71 Rj. The RV detection by Sato 
et al. (2005) was followed by discovery of the 0.3% deep 
transits by the same team. 

Among the ~ 20 Neptune-mass objects found by RV 
searches, only one is known to transit. This is GJ 436b, 
whose transits were recovered by Gillon et al. (2007). 
GJ 436b is thus an extremely valuable and unique ob- 
ject, the only Neptune-mass planet other than our own 
Uranus and Neptune, where the radius has been deter- 
mined (R w 4.9 R e or 4.2 i?®; Bean et al. 2008; Torres, 
Winn & Holman 2008, respectively, hereafter B08 and 
T08), and its internal structure investigated (e.g. Baraffe 
et al. 2008). Based on these results, GJ 436b is a super- 
Neptune with ~ 22 M® total mass with extreme heavy 
element enrichment, and only <~ 10% mass contribution 
by a H/He envelope. 

One of the wide-field surveys involved in the detec- 
tion of TEP's is the HATNet survey (Bakos et al. 2002, 
2004), which currently operates 6 small fully-automated 

13 Of course, the first and most prominent one being HD 209458b 
(Charbonneau et al. 2000; Henry et al. 2000). 



wide-field telescopes. One station is the Fred Lawrence 
Whipple Observatory (FLWO) of the Smithsonian Astro- 
physical Observatory (SAO) on Mt. Hopkins in Arizona 
with four telescopes (HAT-5, HAT-6, HAT- 7, HAT-10), 
and the other is the rooftop of the Submillimeter Array 
Hangar (SMA) of SAO atop Mauna Kea, Hawaii. These 
telescopes are modest 0.11m diameter f/1.8 focal ratio 
tclcphoto lenses that arc using front-illuminated CCDs 
at 5-min integration times. 

Here we report on HATNet's discovery of HAT-P-llb, 
the second transiting hot Neptune known. HAT-P-llb 
orbits the bright (V=9.59) K dwarf star GSC 03561- 
02092 (hereafter called HAT-P-11) with a period of 
4.8878 days, and when it transits the star it causes a dip 
in the stellar light curve of about 4.2millimag (mmag). 
HAT-P-llb is certainly the smallest radius planet found 
by ground-based transit searches, the only planet known 
with smaller radius being Corot-7b (Leger et al. 2009). 
Importantly, the coordinates of the parent star place it 
on one of the detectors of the forthcoming Kepler mis- 
sion; this should allow a broad range of useful follow-on 
observations to characterize both the planetary system 
and the parent star. Because HAT-P-11 is a bright star, 
but still conveniently below the bright limit of the Kepler 
mission, the extraordinary precision of repeated mea- 
surements made over several years should lead to very 
accurate characterization of the system. We note that 
HAT-P-llb would have not been detected by HATNet if 
it orbited a significantly earlier star, such as the typical 
F and G dwarf stars making up the bulk of the HATNet 
transit candidates. 

The layout of the paper is as follows. First we de- 
scribe in § 2 the observational data that led to the discov- 
ery of HAT-P-llb, including the photometric discovery 
data, the reconnaissance spectroscopic observations, the 
photometric follow-up, and the high resolution and high 
signal-to-noise (S/N) spectroscopy. Then we determine 
the parameters of the host star HAT-P-11 (GSC 03561- 
02092) in § 3 by exploring a number of alternate ways. In 
§ 4 we investigate whether the observational data is due 
to a system that mimics planetary transits, and prove 
that this is not the situation, and HAT-P-llb is a bona- 
fide planet. We go on in § 5 to perform global model- 
ing of the data to determine system parameters, such 
as the orbital and planetary parameters. We elaborate 
on treating the systematic variations in an optimal way, 
and present a full analysis that takes them into account 
(§ 5.2, § 5.1 and the Appendix). Finally, we discuss the 
implications of our findings in the Discussion (§ 6). 

2. OBSERVATIONS 

2.1. Photometric Discovery 

The 8.4° x 8.4° region around GSC 03561-02092, a field 
internally labeled as "G155", was observed on a nightly 
basis in two seasons, whenever weather conditions per- 
mitted. First, during the Fall of 2004, we acquired 1213 
frames with the HAT-6 instrument located at FLWO, 
and 4091 frames with HAT-9 located at SMA, Mauna 
Kea. We revisited the field in 2005, and acquired an addi- 
tional 6166 frames with HAT-9. Altogether we gathered 
11470 5-min exposures at a 5.5-minute cadence. This 
unusually rich dataset was motivated by the overlap of 
G155 with the field of view of the future Kepler mission. 

The calibration of the HATNet frames was done uti- 



HAT-P-llb 



3 



8.345 r j 

8 8.350 '- 

| 8.355 - - 
E : 

I 8 360 ' U : 

o) 8.365 - - 
5 8.370 '- 

8.3/5 r " 

-0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 
Orbital phase 

Fig. 1.— The HATNet discovery light curve of HAT-P-11 ex- 
hibiting 11470 individual measurements at 5.5 min cadence. The 
unbinncd instrumental /-band photometry was obtained with the 
HAT-6 (Arizona) and HAT-9 (Hawaii) telescopes of HATNet (see 
text for details), and folded with the period of P = 4.8878162 days, 
which is the result of the global fit described in § 5. Zero orbital 
phase corresponds to the center of the transit. Superimposed is 
the so-called "P1P3" analytic model that was used to describe the 
HATNet data (an approximation of the Mandcl & Agol (2002) an- 
alytic formulae; see § 5.1). 

lizing standard procedures based on IRAF . The cali- 
brated frames were then subjected to star detection and 
astrometry, as described by Pal & Bakos (2006). Aper- 
ture photometry using three apertures was performed 
on each image at the fixed positions of the stellar cen- 
troids, as derived from the 2MASS catalog (Skrutskie et 
al. 2006) and the individual astromctric solutions relating 
the pixel coordinates to the world coordinate (ICRS) sys- 
tem of 2MASS. We extracted photometry for all 125,000 
stars down to / < 14 that fell in the field. The raw in- 
strumental magnitudes m r of each individual frame were 
transformed to a reference frame by fourth order poly- 
nomials in X and Y, and first order in color. The fitted 
magnitudes m / yielded by the above smooth fit are used 
for generating the time versus mj light curves. 

These light curves have a noise characteristic that is 
sometimes referred to as pink-noise (Pont, Zucker, & 
Queloz 2006), because it is a combination of Gaussian 
white noise (due to e.g. photon noise) and a correlated 
red-noise (due to e.g. spatial drift of stars and uncor- 
rected flatfield effects) . The red-noise is often referred to 
as trends or systematic variations. Trends in the light 
curves have an adverse effect on detecting shallow tran- 
siting signals: they mask the real signal, and when phase- 
folded with a trial period they can also mimic transit sig- 
nals. Furthermore, for photometric follow-up data, the 
trends can lead to large systematic errors in the transit 
parameter determination, such as in the impact parame- 
ter or the depth of the transit. Thus, proper treatment of 
systematic variations is crucial for both discovery and ac- 
curate follow-up characterization, especially for the shal- 
low (4.2 mmag) transit of a hot Neptune presented in 
this paper. The two basic methods we have employed 
are the External Parameter Decorrelation (EPD) tech- 
nique, briefly described earlier in Bakos et al. (2007b), 
and the Trend Filtering Algorithm (TFA; Kovacs et al. 
2005). The technical details and some new definitions 
are given in the Appendix (§ A). 

14 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



The HATNet light curves were decorrelated against 
trends using the EPD technique (in constant EPD mode) , 
and subsequently by a simple global TFA (without recon- 
struction, and separately from the EPD). For the bright- 
est stars in the field we achieved a photometric preci- 
sion of 2.9 mmag at 5.5 min cadence. The light curves 
were searched for periodic box-like signals using the Box 
Least Squares transit-search method (BLS, see Kovacs et 
al. 2002). The BLS frequency spectrum of GSC 03561- 
02092 (also known as 2MASS 19505021+4804508; a = 
19 h 50 m 50.21 s , <5 = +48 d 04 m 50.8 s ; J2000) showed a num- 
ber of significant peaks, the highest one at 0.03148 [c/d] 
(31.76d period), and the second one at ~ 0.2046 [c/d] 
(4.8878 days). Fortunately, we inspected the second 
peak, and found the corresponding transit (Fig. 1) with 
a depth of 4.2 mmag, worthy of follow-up. The dip 
had a relative duration (first to last contact) of q w 
0.0206, equivalent to a total duration of Pq w 2.41 hours. 
The signal was confirmed with subsequent reconstruc- 
tive TFA using a trapeze-shaped model function, be- 
ing the most significant periodicity in the data (the 
31.76 d signal was probably due to a left-over system- 
atic). The rms of the TFA processed light curve was 
<~ 3.2 mmag, making the detection of the 4.2 mmag tran- 
sit dip rather challenging, but still robust with S/N = 
11.4. To our best knowledge, this is the shallowest tran- 
sit signal found by ground-based transit searches that 
belongs to a confirmed planet. HD 149026, with an even 
shallower transit signal (Winn et al. 2008) was first de- 
tected via spectroscopy (Sato et al. 2005). The second 
shallowest planetary transit found by the transit search 
method is HAT-P-2b (Bakos et al. 2007b). We note that 
on occasion HATNet has found even shallower transiting 
signals (down to 2 mmag depth) where the transit was 
real, but was due to a blended system. 

2.2. Follow-up reconnaissance spectroscopy 

Following the procedures described in Latham et al. 
(2008), we used the CfA Digital Speedometers (DS; 
Latham 1992), mounted on the 1.5-m Wyeth Reflector at 
the Oak Ridge Observatory in Harvard, Massachusetts, 
and on the 1.5-m Tillinghast Reflector at FLWO, to ob- 
tain low S/N ratio high-resolution spectra. Altogether 5 
spectra were obtained between 2001 December 17 and 
2007 June 5 (the first observations were obtained for 
a different project to survey dwarfs in the solar neigh- 
borhood that might be suitable targets for SETI). The 
signal-to-noise ratios ranged from 10 to 18 per spectral 
resolution element of 8.5 km s -1 . 

Reconnaissance spectroscopy is an important step in 
weeding out astrophysical false positive systems that 
mimic planetary transits. Stellar parameter determina- 
tion can distinguish between dwarf and giant stars via 
the measurement of surface gravity logg*, and thus elim- 
inate systems where the light curve of an eclipsing binary 
is blended with a giant star, or where the transit signal 
seen in our data may not be real, as it is physically not 
feasible to orbit around (and outside) the giant star with 
such a short period. In addition, large RV variations 
(of the order of several kms -1 ) are indicative of orbital 
motion due to stellar companions rather than planets or- 
biting a star. Fine analysis of the spectra can also reveal 
stellar triple systems. Finally, the rapid rotation of a host 
star, as may be indicated by the rotational broadening of 
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TABLE 1 

Photometric follow-up observations of 
HAT-P-11 



Local Date 


Instrument 


Filter 


Ntr 


Type 


2007.0902 


FLW012 


z 





0IBE0 


2007.0907 


FLW012 


z 


1 


--BE0 


2007.1021 


FLW012 


z 


10 


0IBE(0) 


2008.0518 


FLW012 


z 


53 


(OI)BEO 


2008.0603* 


Schmidt 


Ic 


56 


(DIB--) 


2008.0701 


FLW012 


z 


62 


OIB— 


2008.0829* 


Schmidt 


Ic 


74 


(0IBE0) 


2008.0903 


Schmidt 


Ic 


75 


0IBE0 


2008.1002 


FLW012 


z 


81 


OIBEO 


2008.1007 


FLW012 


z 


82 


OIBEO 


2008.1115 


FLW012 


z 


90 


OIB — 


2008.1120 


FLW012 


r 


91 


OIBEO 



Note. — Dates marked with "*" were of poor 
quality, and not used in the analysis. Ntr gives the 
transit event number, counted from the 0th follow-up 
event on 2007 September 2. The Type column shows 
which parts of the transit were caught: Out-of-transit 
(OOT), Ingress, Bottom, Egress, OOT. Parentheses 
mark marginal data quality. 

the spectra, is often correlated with a stellar companion 
that is massive enough to synchronize the rotation with 
the orbital motion. 

HAT-P-11 survived all these steps, and the RV mea- 
surements showed an rms residual of 0.29 km s -1 , con- 
sistent with no detectable RV variation. Initial at- 
mospheric parameters for the star, including the effec- 
tive temperature T^ff* = 4750 ± 125 K, surface gravity 
log = 4.5 ± 0.25 (cgs), and projected rotational ve- 
locity vsini — 0.0 km s -1 , were derived as described by 
Latham et al. (2008). The mean line-of-sight velocity of 
the star was T = — 63.56 ± 0.29 kms -1 on an absolute 
scale. 

2.3. Photometric follow-up observations 

Photometric follow-up observations are important i) 
to perform independent confirmation of the initial detec- 
tion, ii) to allow for accurate characterization of the sys- 
tem (see also § 5), iii) to help eliminate blend scenarios 
(§ 4), and, if multiple events are observed, iv) to possi- 
bly average out extra variations of the light curve due 
to spots. Since the transit of HAT-P-11 was originally 
detected with the 11cm diameter HATNet telescopes (al- 
though when phase-folding many events), it is certainly 
feasible to confirm a single transit with a lm-class tele- 
scope. However, acquiring a high-quality observation of 
the 4.2mmag transit dip in the light of a star as bright 
as V = 9.59 is quite challenging. An accurate and highly 
precise light curve is essential for eliminating blend sce- 
narios, and then for determining the physical parameters 
of the transiting planet-star system. In general, the shal- 
lower the transit, the wider the range of possible blends 
(i.e. more than two body systems) that can mimic the ob- 
served transit. Some of these can only be distinguished 
by subtle effects, such as the duration of ingress/egress, 
and accurate depth and shape of the transit. 

For these reasons, we launched an extensive photome- 
try follow-up campaign, and attempted transit observa- 
tions of HAT-P-llb altogether 12 times in 2007 and 2008, 
leading to ~ 10 successful observations of partial or full 
transits (see Tab. 1). We primarily used the FLWO 1.2 m 
telescope and the KeplerCam CCD in Sloan z-band, with 
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Fig. 2. — Unbinncd instrumental transit photometry follow-up 
light curves acquired by the KeplerCam at the FLWO 1.2 m tele- 
scope (in Sloan z band, and Sloan r band) and by the 0.6 m Schmidt 
telescope at Piszkesteto, Hungary (in Iq band). Superimposed is 
the best-fit transit light curve model, yielded by the global analy- 
sis described in § 5. The light curves have been subject to global 
EPD and TFA corrections using an analytical transit model for the 
intrinsic signal (sec § A). 

exposure times of ~ 10 sec and read-out time of 12 sec. 
The last observation on 2008 November 20/21 MST was 
taken through Sloan r band to complement our blend 
analysis (§ 4). We also observed transits 3 times using 
the 0.6 m Schmidt telescope of Konkoly Observatory at 
the Piszkesteto mountain station. These observations 
were taken through J^-band, and the only useful dataset 
proved to be the one from 2008 September 3/4 CET. 

Data were reduced in a similar manner as for the 
HATNet data (§ 2.1), and as described in Bakos et al. 
(2008). Following bias and flat calibration, we derived 
an initial first order astrometrical transformation be- 
tween the ~ 750 brightest stars and the 2MASS cata- 
log. The position of HAT-P-11 in the 2MASS catalog 
was corrected for each observed epoch before deriving 
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the astromctry, due to its moderately high proper motion 
(263. 32± 1.31masyr _1 ; Perryman et al. 1997). Photom- 
etry was carried out for all stars in the field using three 
apertures that were adjusted each night to match the ob- 
serving conditions (sky background, profile width). For 
the nights with good seeing, one aperture was kept small 
to avoid the faint neighbor 2MASS 19505049+4805017 
currently at approx 6" distance (see § 4). Instrumental 
magnitudes were transformed to a photometric reference 
frame (selected to be at low airmass, low sky background, 
etc.) using a first order polynomial in the X, Y pixel co- 
ordinates and the 2MASS J-K color of ~ 600 stars. The 
transformation was iteratively determined by weighting 
with individual Poisson noise errors of the stars, and re- 
jecting 3-cr outliers, plus eliminating the main target and 
its faint neighbor from the fit. The smooth fit was re- 
peated based on the ~ 180 best stars, using the median 
of the individual light curve magnitude values as a new 
reference system, and weighting the fit with the rms of 
the light curves (i.e. substituting the former Poisson er- 
rors). 

The FLWO 1.2 m KeplerCam observations usually re- 
sult in high quality photometry of ~1% deep planetary 
transits (see e.g. Latham et al. 2008; Bakos et al. 2008), 
because of the large field-of-view (FOV = 23' x 23') 
with many potential comparison stars, fine pixel res- 
olution (0.336"pixcl~ 1 ), good quality and high quan- 
tum efficiency CCD (monolithic, 4K x 4K Fairchild 486 
chip), good sky conditions from FLWO, the fast readout 
through 4 channels, and careful data processing. The 
performance on this very shallow transit of a fairly bright 
star, however, turned out to be slightly sub-optimal with 
some residual trends, necessitating diverse application of 
the EPD and TFA methods. Since these were part of our 
global modeling of the data, including the HATNet and 
RV data, they are detailed later in § 5 and in the Ap- 
pendix. The follow-up light curves after applying simul- 
taneous EPD-TFA with per-night coefficients for EPD 
and global coefficients for TFA are displayed on Fig. 2. 
The model function (denoted as mo(p,ti) in Eq. Al in 
the Appendix) in the simultaneous EPD-TFA fit is an 
analytical transit model from Mandcl & Agol (2002). 

2.4. High resolution, high S/N spectroscopy 

We started observations of HAT-P-11 on 2007 August 
22 with the HIRES instrument (Vogt et al. 1994) on the 
Keck-I telescope located on Mauna Kea, Hawaii. It was 
soon realized that the target is unusually complicated 
because of the small RV amplitude with respect to the 
moderate velocity jitter due to the active K dwarf star, 
and also because of long term trend(s) present. Thus, 
HAT-P-11 has been extensively observed over the past 
two years, and we have gathered altogether 50 spectra 
and 3 template observations (see Tab. 3). This is 5- 
10 times the number of spectra collected for a typical 
HATNet transit candidate. 

The width of the spectrometer slit used on HIRES was 
CC'86, resulting in a resolving power of A/AA w 55,000, 
with a wavelength coverage of ~ 3800 — 8000 A. The io- 
dine gas absorption cell was used to superimpose a dense 
forest of I2 lines on the stellar spectrum and establish an 
accurate wavelength fiducial (see Marcy & Butler 1992). 
Relative RVs in the Solar System barycentric frame were 
derived as described by Butler et al. (1996), incorporat- 



ing full modeling of the spatial and temporal variations 
of the instrumental profile. The final RV data and their 
errors as a function of barycentric Julian date are listed 
in Tab. 3. It is reassuring that a simple Fourier analysis 
of the RV data without prior assumption on any periodic 
event yields a primary period of 4.888 ± 0.019 days, con- 
firming the same periodic phenomenon that is present in 
the independent dataset consisting of the discovery and 
follow-up photometry. 15 

Based on the numerous observations of a transit-like 
event in the photometry, our initial physical model was 
thus a single body orbiting and transiting the star, caus- 
ing both the RV variations and the transits. The RV 
variations induced on the parent star can be character- 
ized by six parameters: the period P, the center of transit 
T c (i.e. the phase of the RV curve), the semi-amplitude 
K , the RV zero-point 7, and the Lagrangian orbital el- 
ements (k,h) — e x (cosu, sinw). An orbital fit to the 
RV data without any constraints from the photometry 
yielded the following values: P = 4.8896 ± 0.0017 d, 
T c = 2454552.04 ± 0.16 (BJD), K = 13.1 ± 2.1 ms" 1 , 
7 = 0.9±0.8ms _1 (arbitrary scale), k = 0.28±0.12 and 
h = 0.09 ±0.09. 

Later (§ 5) we present a global modeling of the data, 
where joint analysis of the photometry and RV data is 
performed. This yields precise P and T c parameters that 
are almost entirely constrained by the photometric data 
(notably the sharp ingress/egress features), with virtu- 
ally no coupling from the RV data. The knowledge of 
period and phase from photometry is a very tight con- 
straint in the analysis of the RV data, yielding signifi- 
cant improvement on determination of the other orbital 
parameters. It is thus well justified to analyze the RV 
data separately by fixing the ephemeris to that deter- 
mined by the photometry, and assuming that there is an 
eccentric orbital motion present. In addition, we can also 
investigate whether our initial model correctly described 
the physics of the system, and whether there are other 
signals present in the RV data. 

Along these lines, our refined model for the RV data 
was an eccentric orbital motion with P and T c fixed to 
those values found from the global modeling (and pri- 
marily constrained by the photometry) , plus a sinusoidal 
motion of the form A 2 sm(f 2 t) + B 2 sin(/2i) (equivalent 
to A' 2 sin(/ 2 i + 0), but linear in the fitted A 2 and B 2 pa- 
rameters). We searched the fi domain by fitting 7, K, k, 
h (eccentric orbit parameters) and A 2 , B 2 at each f 2 , and 
noting the \ 2 . This "frequency scan" located a number 
of significant peaks with small inverse x 2 (Fig- 3, upper 
panel), the most significant being f 2 — 0.9602 ± 0.0009. 
We suspected that the peaks around 1 [c/d] might be due 
to a long-term trend in the RV, sampled with a daily and 
lunar-cycle periodicity. To test this effect we generated 
a mock signal by co-adding an analytic eccentric orbital 
motion, a long term RV drift, and Gaussian noise, and 
sampled these at the exact times of the 50 Keck observa- 
tions. We ran the same fitting procedure (Keplerian plus 
sinusoidal motion) in the above described frequency scan 
mode on the mock data. Indeed, the drift appeared as a 
number of significant frequencies around 1 [c/d] (Fig. 3, 
panel b), their exact location depending on the ampli- 

15 The uncertainty of the period was simply derived from the full 
width at half magnitude of the appropriate peak in the spectrum. 
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TABLE 2 



Follow-up Photometry for HAT-P-11 


BJD 


Mag 


Error 


Raw Mag 


Filter 


Mr 


Instrument 


c~i AC A1 A C C-\ T A C 

2454o4o.ol74o 


-0.00134 


0.000b4 


-0.00082 


z 





T7T Tli 7/— \ 1 O 

r LWU12 


2454o4o. 61771 


-0.00228 


n /"inn/" 1 i 

0.000b4 


-0.00234 


z 





r LWU12 


2454351. 74387 


0.0U545 


O.OOObb 


0.00bl5 


z 


1 


T7T 7/— * 1 O 

r LWU12 


2454ool. 74414 


0.00320 


O.OOObb 


0.00134 


z 


1 


T7T 7/— * 1 O 

r LWU12 


2454o95. 58419 


-0.00509 


0.000b8 


-0.00078 


z 


10 


TTT ^ Jl 7/— * 1 O 

r LWU12 


2454395.58575 


-0.00042 


0.000b2 


-0.00170 


z 


10 


7/— \ 1 c 1 

I 1 LW012 


2454oUo.o49U5 


0.00287 


0.000b5 


0.00393 


z 


53 


TT1T TT 7/— 1 1 O 

r LWU12 


2454b05. 84933 


0.00280 


n /"inn/" 1 rr 

0.000b5 


0.007b3 


z 


53 


T7T 7/— * 1 O 

r LWU12 


2454b49. 8171b 


0.00186 


n r\ c\ r\ a c\ 

0.000b9 


0.00104 


z 


62 


r LWU12 


2454649.81742 


0.00168 


0.00009 


-0.00059 


z 


62 


FLW012 


2454713.34579 


-0.00234 


0.0009b 


0.00018 


I 


75 


Schmidt 


2454713.34612 


-0.00136 


0.00099 


0.00131 


I 


75 


Schmidt 


2454742.67630 


-0.00239 


0.00053 


-0.00358 


z 


81 


FLW012 


2454742.67664 


-0.00049 


0.00053 


-0.00093 


z 


81 


FLW012 


2454747.57949 


0.00194 


0.00054 


-0.00019 


z 


82 


FLW012 


2454747.57981 


-0.00251 


0.00054 


-0.00376 


z 


82 


FLW012 


2454786.57057 


0.00386 


0.00001 


0.00543 


z 


90 


FLW012 


2454786.57106 


0.00042 


0.00001 


0.00307 


z 


90 


FLW012 


2454791.55380 


0.00134 


0.00050 


0.00048 


r 


91 


FLW012 


2454791.55562 


-0.00340 


0.00050 


-0.00356 


r 


91 


FLW012 



Note. — column (1): Barycentric Julian Day, 
column (2): Best detrended magnitude, normalized to 0.0 out of transit, 
column (3): Estimated error in the best magnitude, 
column (4): Magnitude before detrending (denoted as nif in the text), 
column (6): Transit number. See Tab. 1. 

This table is presented in its entirety in the electronic edition of the Astro- 
physical Journal. A portion is shown here for guidance regarding its form and 
content. 



tudc of the drift and the random seed used to generate 
the noise. Then, we repeated the above analysis includ- 
ing a long term drift in the fitted model function while 
performing the frequency scanning procedure (i.e. fitting 
at each frequency value that is being stepped on a grid) . 

The results are plotted in the c) panel of Fig. 3, on 
the same scale as used in panel a) and b). It can be 
clearly seen that the magnitude of the residuals is defi- 
nitely smaller and the strong structures disappeared. Us- 
ing the mock data set, the results are almost the same 
(Fig. 3 panel d), confirming our assumptions that a long 
term drift in the data causes aliases around 1 [c/d]. This 
is in line with the experience based on long-term mon- 
itoring of stars with Kcck/HIRES: in nearly every case 
where there is either a planet or a long-term trend, a 
spurious spike appears in the periodogram at close to 1 
day. 

Note, however, that in spite of the clearing of peaks 
around l[c/d], the 0.9602 [c/d] peak remained present 
even after the simultaneous fit of a Kcplerian orbit and 
a drift to the Keck data. 

It is also interesting that the x 2 of the fit with Kep- 
lerian plus f 2 sinusoid component is somewhat smaller 
than the fit with Keplerian plus long term trend. The 
fi = 0.9602 [c/d] periodicity may be a real physical ef- 
fect, or can be an alias of a real physical effect with differ- 
ent period (rotation period, activity of star), or an alias 
of a systematic effect (lunar cycle, instrumental effect). 
Altogether we have three components to consider in the 
final model of the RV variations: i) Kcplerian orbit, ii) 
long-term drift, iii) f 2 = 0.9602 [c/d] periodicity. 

The model of a Keplerian orbit and a long-term drift 
(i+ii) has simple underlying physics, by assuming an in- 
ner and an outer planet. The results of the orbital fit 
for this basic model are exhibited on Fig. 4. The data 




Frequency 



Fig. 3. — (Top panel:) The inverse x 2 as a function of fre- 
quency when fitting a combined Kcplerian orbit with fixed P and 
T c ephemeris, plus a sinusoidal component with f'2 frequency to 
the RV observations. The high peaks indicate good fit with small 
X 2 - The 0.96[c/c(] peak is marked with a circle. (Second panel:) 
The same plot, but on a mock dataset that consists of a Keple- 
rian, a drift, and Gaussian noise. (Third panel:) Same as the top 
panel, but the fitted model was extended by a drift with fitted 
slope. Much of the 1 [c/d] frequencies disappeared, but the one at 
e.g. 0.9602 [c/d] remained. (Bottom panel:) Same as the second 
panel (mock data), with fitted model extended with a linear slope. 
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folded with the P = 4.8878162 d period after removal 
of the best-fit drift is plotted in the top panel. Superim- 
posed is the Keplerian orbit that is clearly eccentric. The 
orbital elements for this fit were: K = 11.6 ± 1.2ms -1 , 
k = ecosw = 0.201 ± 0.049, h = esinw = 0.051 ± 0.092, 
7 = — 0.4 ± 0.8 ms- 1 , and the best fit value for the linear 
drift is Gi = 0.0297± 0.0050ms- 1 day" 1 . Note that the 
K amplitude of the orbit is only 11.6±1.2ms . In order 
to have a reduced x 2 value of 1.0, a jitter of 5.01ms- 1 
has to be added to the formal errors. 

We also derived an orbital fit with a full model of 
all three components (Keplerian, drift, short-period si- 
nusoidal) so as to check the effect of /2 on the or- 
bital parameters. We found that the change in the or- 
bital elements is insignificant, with small decrease in 
their respective error-bars: K = 10.91 ± 0.96ms- 1 , 
k = ecosw = 0.225± 0.038, h = esinw = 0.069± 0.072, 
drift Gi = 0.0176 ± 0.0056ms- 1 day" 1 , and the ampli- 
tude of the /2 sinusoid is 5.7± 1.3ms- 1 . The jitter value 
changed to 3.89 ms" 1 . 

We also investigated models with non- linear drift, as 
characterized by G2 quadratic and G3 cubic terms, and 
found that these are insignificant based on the present 
data. We also checked for correlations between the RV 
residuals from the best fit, and the spectral bisector- 
spans, and the S activity index, but found that these 
correlations were insignificant. 

Given the fact that the origin of the ji periodicity 
(component iii) is unknown, and it is suspected to be 
an alias that may diminish by taking more data, plus 
the key orbital elements do not change significantly by 
taking it into account, in the rest of this paper we 
adopted the simpler model of a Keplerian orbit (with 
K = 11.6 ± 1.2 ms- 1 ) and a long-term drift without the 
short-term periodicity. 

3. PROPERTIES OF THE PARENT STAR 

Knowledge of the parameters of the host star is im- 
portant, because it puts the relative quantities arising 
from the global modeling of the photometric and RV 

data (§ 5) on an absolute scale, since the planetary ra- 

2/3 

dius R p is cx i?* and the planetary mass M p is cx MJ . 
Also, modeling of the photometric transit requires limb- 
darkening parameters. These may be fitted along with 
other parameters if the data are of high quality and al- 
low it (Brown et al. 2001). Alternatively, limb-darkening 
coefficients are taken from look-up tables, such as that 
of Claret (2004), which depend on stellar atmospheric 
parameters, primarily effective temperature. Conversely, 
the presence of the planetary transit imposes constraints 
on the stellar parameters through the normalized semi- 
major axis a/R± (and stellar density p*) as yielded by 
the global modeling of the photometry, and, for eccen- 
tric orbits, the RV data (T08). 

3.1. Basic Stellar Parameters 

We employed the Spectroscopy Made Easy (SME) 
package of Valenti & Piskunov (1996) along with the 
atomic-line database of Valenti & Fischer (2005) to de- 
rive an initial value for the stellar atmospheric param- 
eters. In this analysis we used the iodine-free template 
spectrum obtained by the HIRES instrument on Keck I. 
The SME analysis of stellar spectra resulted in a stel- 
lar surface gravity logg* = 4.7 ± 0.1 (CCS), metal- 



20 




Fig. 4. — (Top panel:) The radial velocity measurements from 
Keck for HAT-P-11 after correction for a linear drift seen in the 
data. The measurements are folded using our best-fit ephemeris 
from global modeling of the photometry and RV data (see § 5), 
and superimposed by our (eccentric) orbital fit. The center-of- 
mass velocity has been subtracted. The orbital phase is shifted 
to be zero at the center of the transit. The error-bars have not 
been inflated with the jitter (5.01ms _1 ). The dashed line (barely 
discernible from the thick line) is the orbital fit with the same 
ephemeris and K semi-amplitude, but using the refined kc and he 
Lagrangian orbital elements (§ 5.5). (Second panel:) RV residuals 
after the orbital fit. (Third panel:) The S activity index of HAT-P- 
11 phase-folded. (Bottom panel:) Bisector spans (BS) for the Keck 
spectra plus the three template spectra, computed as described in 
the text. The mean value has been subtracted. The vertical scales 
on the lower and upper panels are the same. 

licity [Fe/H] = 0.32 ± 0.06 dex, effective temperature 
T e s* — 4850 ± 50 K, and the projected rotational ve- 
locity vsmi — 0.5 ± O.Skms- 1 . However, the vsini 
value depends inversely on the assumed value of macro- 
turbulcncc v mac in the star. The value of v mac assumed 
in the SME derivation was 2.57kms _1 , but for stars of 
type K0 or later, v mac may be as small as 1.3kms _1 
(e.g. Gray 1988). With that value of v mac , vsini could 
be as high as 2.7kms- 1 . Because of the indeterminacy 
of v mac for HAT-P-11, we can conclude solely that HAT- 
P-11 is a relatively slow rotator, and assign it a projected 
rotational velocity vsini = 1.5 ± 1.5 km s -1 , which en- 
compasses both extremes of usinz noted here. 

At this stage we could use the effective temperature 
as color indicator and the surface gravity as luminosity 
indicator, and determine the stellar parameters based on 
these two constraints using a set of isochrones. However, 
it has been shown (Sozzetti et al. 2007) that log g+ has 
a subtle effect on the spectral line shapes, and is usually 
not the best luminosity indicator. For planetary tran- 
sits, the a/R+ normalized semi-major axis and related 
p+ mean stellar density typically impose a stronger con- 
straint on possible stellar models (Sozzetti et al. 2007). 
(The use of a/R* is discussed later in § 3.2.) However, 
for HAT-P-11 there is an even better luminosity indi- 
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TABLE 3 

Relative radial velocity measurements of HAT-P-11 



BJD RV cr RV O-C BS o-bs S 

(ms _1 ) (ms _1 ) (ms _1 ) (ms _1 ) (ms _1 ) 



2454335.89332 


8.31 


1 


10 


1.82 


-2.40 


5.32 


0.603 


2454335.89997 


8.61 


1 


.17 


2.03 


-2.70 


5.38 


0.604 


2454336.74875 


-3.12 


1 


.12 


-6.29 


-17.80 


6.29 


0.614 


2454336.25715 










-1.30 


5.27 


0.616 


2454336.86162 


-0.34 


1 


11 


-1.36 


-1.74 


5.41 


0.613 


2454336.94961 


-1.60 


1 


13 


-0.93 


19.68 


6.98 


0.621 


2454337.73150 


-8.74 


1 


.15 


2.78 


21.40 


7.28 


0.641 


2454337.92101 


-15.25 


1 


.11 


-2.42 


8.67 


5.78 


0.628 


2454338.74219 


-14.18 


1 


.10 


0.02 


25.38 


7.47 


0.645 


2454338.92046 


-15.82 


1 


.06 


-2.09 


18.73 


6.63 


0.622 


2454339.89292 


-6.53 


1 


10 


-0.09 


-5.27 


5.60 


0.637 


2454343.83879 


-11.99 


1 


10 


1.48 


1.34 


5.47 


0.623 


2454344.96843 


— 11.31 


1 


.07 


-7.44 


2.34 


5.44 


0.610 


2454344.47643 










-9.38 


5.69 


0.616 


2454344.48214 










-9.88 


5.57 


0.612 


2454396.75318 


-1.40 


1 


.00 


10.52 


-5.31 


5.69 


0.641 


2454397.76902 


-5.30 


1 


,04 


5.85 


4.01 


5.51 


0.638 


2454427.72779 


2.75 


1 


.26 


8.18 


7.61 


5.87 


0.592 


2454428.69600 


12.16 


1 


.39 


3.77 


-6.33 


5.88 


0.580 


2454430.69738 


-4.82 


1 


22 


4.64 


-4.12 


5.52 


0.578 


2454454.70749 


-4.02 


1 


.58 


0.48 


-9.30 


6.36 


0.564 


2454455.69785 


-8.06 


1 


.39 


2.78 


-10.02 


6.20 


0.581 


2454456.68976 


-6.75 


1 


.47 


1.02 


1.94 


5.72 


0.559 


2454457.68670 


3.02 


1 


.39 


-0.98 


-8.06 


5.55 


0.568 


2454460.72727 


-5.84 


1 


.51 


4.91 


0.15 


5.72 


0.611 


2454549.09456 


-18.39 


1 


.66 


-10.98 


9.97 


5.83 


0.540 


2454549.13273 


-21.19 


1 


.65 


-13.92 


11.64 


5.93 


0.541 


2454602.01197 


-13.45 


1 


10 


-7.82 


2.08 


5.42 


0.591 


2454602.11819 


-12.58 


1 


.31 


-6.57 


-2.61 


5.25 


0.585 


2454602.94861 


-4.49 


1 


19 


0.97 


3.85 


5.46 


0.591 


2454603.00783 


-6.92 


1 


.23 


-1.71 


3.64 


5.54 


0.593 


2454603.11953 


-6.07 


1 


16 


-1.42 


-5.07 


5.78 


0.588 


2454603.95329 


7.05 


1 


.34 


4.05 


1.34 


5.40 


0.589 


2454604.03295 


3.47 


1 


.35 


-0.60 


1.08 


5.37 


0.595 


2454604.12769 


3.48 


1 


.21 


-1.95 


6.46 


5.65 


0.595 


2454634.00592 


14.22 


1 


.27 


-0.66 


-12.14 


6.68 


0.582 


2454634.95641 


18.14 


1 


.14 


6.98 


-2.23 


5.52 


0.596 


2454636.01011 


1.45 


1 


.16 


4.90 


-0.31 


5.39 


0.598 


2454637.99531 


2.11 


1 


.13 


0.21 


2.64 


5.55 


0.623 


2454639.02952 


15.89 


1 


.08 


-0.81 


-1.66 


5.79 


0.613 


2454640.07995 


2.31 


1 


.22 


-4.53 


-1.84 


5.73 


0.611 


2454641.05878 


-8.01 


1 


.24 


-3.80 


3.15 


5.49 


0.614 


2454674.98972 


-9.75 


1 


.33 


-8.38 


-3.81 


5.67 


0.587 


2454722.80946 


18.87 


1 


.16 


2.70 


5.44 


5.86 


0.671 


2454723.80199 


2.07 


1 


.21 


1.40 


-16.64 


6.54 


0.657 


2454725.81649 


2.56 


1 


.21 


-0.29 


5.75 


5.67 


0.662 


2454726.85018 


13.91 


1 


.17 


-3.43 


-10.07 


5.78 


0.657 


2454727.77201 


12.20 


1 


.20 


-2.75 


0.29 


5.57 


0.648 


2454727.90334 


13.77 


1 


.16 


1.32 


-3.79 


5.63 


0.646 


2454777.79897 


2.23 


1 


.21 


1.68 


1.74 


5.59 


0.637 


2454778.78504 


10.39 


1 


.61 


11.00 


-2.29 


5.80 


0.607 


2454779.80416 


12.24 


1 


.33 


5.34 


-5.81 


6.20 


0.632 


2454790.68835 


22.38 


1 


.31 


0.04 


-8.66 


5.76 


0.621 



cator, since it is a bright and nearby star with known 
parallax of small relative error. HAT-P-11 enters the 
Hipparcos catalogue as HIP 97657 with reported paral- 
lax of 27.50 ± 0.96 mas (Perryman et al. 1997), equiva- 
lent to a distance of 36.4 ± 1.3 pc and distance modulus 
of A — 2.80 ± 0.08. Combination of the distance infor- 
mation with the apparent brightnesses in various photo- 
metric bands yields the absolute magnitude of the star, 
which is a tighter constraint on the luminosity than the 
logg* or a/R± constraints. 

As regards apparent magnitude, the TASS (Droege et 
al. 2006) photometry for this star is Vtass = 9.587 ± 
0.071 and J T ass = 8.357 ±0.050, yielding {V - I) T ass = 



1.230 ± 0.087. On the other hand, the Yonsei-Yale stel- 
lar evolution models (Yi et al. 2001; Demarque et al. 
2004) indicate (V - I)YY,cx P cct = 0.98 ± 0.03 for a star 
with T eff * ~ 4850 K (based on the SME results). The 3-cr 
inconsistency cannot be explained by interstellar redden- 
ing, since this is a close-by star. Thus we opted not to 
use the TASS photometry in this analysis. 

The 2MASS catalogue (Skrutskie et al. 2006) pro- 
vided much better agreement. The magnitudes reported 
in the 2MASS catalogue have to be converted to the 
standard ESO system, in which the stellar evolution 
models (specifically the YY models) specify the colors. 
The reported magnitudes for this star are J2MASS = 
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7.608 ± 0.029, H 2 mass = 7.131 ± 0.021 and K 2M ass = 
7.009 ± 0.020; which is equivalent to J = 7.686 ± 0.033, 
H = 7.143 ± 0.028 and K = 7.051 ± 0.021 in the 
ESO photometric system (see Carpenter 2001). Thus, 
the converted 2MASS magnitudes yield a color of (J — 
K) = 0.635 ± 0.043 that is within 1-a of the expected 

(J - JQvY.expect = °- 59 ± °- 02 - Wc thus relicd on tnc 

2MASS K apparent magnitude and the parallax to derive 
an absolute magnitude of M K = 4.18 ± 0.07. The choice 
of K band was motivated by the longest wavelength with 
smallest expected discrepancies due to molecular lines in 
the spectrum of this K4 dwarf. 

In practice, the isochrone search for the best fit stellar 
parameters was done in a Monte Carlo way, by assuming 
Gaussian uncertainties for the Hipparcos parallax, T e ff*, 
[Fe/H], the apparent 2MASS magnitudes, and the con- 
version coefficients by Carpenter (2001) that transform 
the 2MASS magnitudes to the standard K band. Lack- 
ing information, possible correlations between some of 
these parameters (e.g. T e g* and [Fe/H]) were ignored. A 
large set (~ 5000) of random A (distance modulus), T e ff*, 
[Fc/H] and K values were generated, and for each com- 
bination we searched the stellar evolutionary tracks of 
the Yonsei-Yale models for the best fit stellar model pa- 
rameters (such as M*, i?*, etc). For an uncvolved K star 
there was no ambiguity in the solution, i.e. we did not 
enter a regime where isochrones cross each other. Cer- 
tain parameter combinations in the Monte Carlo search 
did not match any isochrone. In such cases (~ 14% of 
all trials) we skipped to the next randomly drawn pa- 
rameter set. At the end we derived the mean values and 
uncertainties of the physical parameters based on their a 
posteriori distribution. We also refined the stellar surface 
gravity. The new value log<?* = 4.59 ± 0.03 agrees well 
with the earlier SME value, and has much smaller un- 
certainty confirming our previous assumptions that the 
absolute magnitude for this star is a better luminosity 
indicator than the surface gravity. 

We then repeated the SME analysis by fixing log g± 
to the new value, and only adjusting T cf f*, [Fe/H] and 
vsini. This second iteration yielded T e g* = 4780 ± 
50 K, and [Fc/H] = +0.31 ± 0.05. It also yielded 
vsini = 0.3 ± 0.5kms _1 , but for the same reason 
as given earlier we adopt the relaxed range v sin i — 
1.5 ± 1.5 km s -1 . We accepted the above values as the 
final atmospheric parameters for this star. We then 
also repeated the isochrone search for stellar parameters, 
yielding M*=0.809;°;o2?M Q , i?*=0.752 ± 0.021 R Q and 
L^=0.26±0.02L Q . Along with other stellar parameters, 
these are summarized in Tab. 4. The stellar evolution- 
ary isochrones for metallicity [Fe/H] =+0.31 are plotted 
in the right panel of Fig. 5, with the final choice of effec- 
tive temperature T c ff* and the absolute magnitude Mk 
marked, and encircled by the 1-a and 2-cr confidence el- 
lipsoids. 

The justification for our choice of using the parallax 
as the luminosity indicator is demonstrated by Fig. 5. 
On the left panel wc plot a set of [Fc/H] =+0.31 Yale 
isochrones as a function of effective temperature, with 
the vertical axis being log <?*, and the observed values 
with their respective error ellipsoids overlaid. The mid- 
dle panel shows the same set of isochrones with the a/i?* 
luminosity indicator on the vertical scale. Here the geo- 
metric semi-major axis a/i?* is determined from the pho- 



TABLE 4 
Stellar parameters for HAT-P-11 



Parameter Value Source 



T eSir (K) 4780 ± 50 SME a 

[Fc/H] +0.31 ± 0.05 SME 

vsini (bus- 1 ) 1.5 ± 1.5 SME 

M* (M ) 0.81±°'° 2 Y 2 +Hip+SME b 

Ri, (Rq) 0.75 ±0.02 Y 2 +Hip+SME 

logs* (cgs) ... 4.59 ±0.03 Y 2 +Hip+SME 

L* (1,0 ) 0.26 ±0.02 Y 2 +Hip+SME 

M F (mag).... 6.57 ±0.09 Y 2 +Hip+SME 

Age (Gyr) .... 6.5^ ? Y 2 +Hip+SME 

Distance (pc) . 38.0 ±1.3 Y 2 +Hip+SME C 



a SME = 'Spectroscopy Made Easy' package 
for analysis of high-resolution spectra Valenti & 
Piskunov (1996). These parameters depend pri- 
marily on SME, with basically no dependence on 
the iterative analysis based on the Hipparcos par- 
allax and YY isochrones (§ 3.1). 

b Y 2 +Hip+SME = Yale-Yonsei isochrones (Yi ct 
al. 2001), Hipparcos distance data, and SME re- 
sults. 

c The distance given in the table is based in the 
self-consistent analysis that relies on the Hipparcos 
parallax and the YY isochrones. It slightly differs 
from the Hipparcos-based distance. 

tometric transit and RV data (as shown later in § 3.2), 
and its error is significantly increased by the uncertain- 
ties in the eccentricity of the RV data. The right panel 
shows the same set of isochrones, with vertical axis being 
the absolute K magnitude. Overlaid is the observational 
constraint based on the apparent magnitude and Hip- 
parcos parallax, along with the respective 1-a and 2-a 
error ellipsoids. It is clearly seen that the right panel 
imposes the tightest constraint on the stellar parame- 
ters. We note here that an isochrone search using the 
reported J — K color (instead of T e ff* as color indicator) 
also agrees with the evolutionary models; however, the 
relative volume of its confidence ellipsoid is somewhat 
larger. 

The effective temperature from the SME analysis and 
the surface gravity derived above correspond to a K4V 
star using Gray (1992). The B — V color index from 
the YY isochrones is 1.063 ± 0.024, also consistent with 
the K4 spectral type of Gray (1992). We also ran the 
Ramirez & Melendez (2005) and Casagrande et al. (2006) 
temperature calibrations in reverse to back out the B — V 
required to produce the SME temperature. We got B — 
V = 1.025 ± 0.023 and 1.067+ 0.025, respectively. For 
final value we accepted their average: B — V = 1.046 ± 
0.024. 

3.1.1. Baraffe isochrones 

We also investigated the dependence of stellar param- 
eters on the choice of isochrones. Since HAT-P-11 is a K 
dwarf, we used the Baraffe et al. (1998) isochrones that 
are usually a better choice for late type dwarfs. Baraffe 
et al. (1998) presented three sets of isochrones with dif- 
ferent mixing length parameters. The one with a = 1.0 
is a better match for low-mass stars (such as late K or 
M dwarfs), the one with a = 1.9 matches the Sun, and 
the third set with a = 1.5 is in between. In Fig. 6 wc 
plot these isochrones as a function of T e fj*, with verti- 
cal axis being absolute K magnitude in the Bessel Brett 
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Fig. 5. — The three panels display the same Yonsei-Yale isochrones between 2.0 and 14.0 Gyrs for metal-rich, [Fe/H] = +0.31 stars, 
including masses between 0.65 Mq < M < 1.25 Mq. The horizontal axis is effective temperature for all cases. The three panels show 
different choices of luminosity indicators (on the vertical axes). The observed values for the stellar parameters of HAT-P-11 are marked by 
the large filled circles, and the 1-cr and 2-cr confidence ellipsoids arc also indicated. Left panel: the luminosity indicator (vertical axis) is 
the stellar surface gravity. Middle panel: The luminosity indicator is the geometric semi- major axis a/R*, as derived from the light curve 
modeling (see § 3.2 and § 5). Right panel: here the luminosity indicator is the Mk absolute K magnitude, based on the 2MASS catalogue, 
transformations by Carpenter (2001), and the Hipparcos parallax. The tightest constraint on the isochrones is provided by the effective 
temperature (based on the SME analysis), and the absolute K magnitude photometry (right panel). 

system. Both the 2MASS K s band and the Baraffe CIT 
systems were transformed to K using the relations in 
Carpenter (2001). The Baraffe isochrones for a > 1.0 are 
given only for solar metallicity, while those for a = 1.0 
are available for solar and sub-solar (Baraffe et al. 1997) 
metallicity only, but not for the metal-rich composition 
of HAT-P-11 ([Fe/H] = +0.31±0.05). To investigate the 
effect of metallicity, we assumed the same qualitative be- 
havior (i.e. opposite shift in the absolute magnitude — 
temperature plane) for isochrones with metallicity in- 
creased by +0.3, independent of their mixing lengths. 
In Fig. 6, we over-plot a metal-poor ([Fc/H] = —0.3, 
a = 1.0) isochrone, and conclude that an opposite change 
in metallicity to match that of HAT-P-11 would move the 
a = 1.0, [Fc/H] = +0.3 isochrone away from the observa- 
tional values of K and T a*. Because of the lack of metal- 
rich models from Baraffe, and the lack of the knowl- 
edge of the proper mixing length parameter, we omit any 
quantitative conclusions from these models. Altogether, 
we found a better match with the YY isochrones, and 
thus accepted YY-based stellar parameters as final val- 
ues (Tab. 4). 
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3.2. Constraints on stellar parameters by the 
normalized semi-major axis 

As noted by Sozzetti et al. (2007); Torres, Winn & 
Holman (2008), a possible luminosity indicator for host 
stars of transiting planets is the a/R* quantity, where 
a is the relative semi-major axis, and i?* is the radius 
of the host star. Here a/i?* is in simple relation with 
the mean stellar density, if we assume that the planetary 
mass is much smaller than the stellar mass: 



Fig. 6. — Baraffe et al. (1998) isochrones with mixing length 
parameters a = 1.0, 1.5 and 1.9 (solid lines) for solar metallicity 
(note that the metallicity of HAT-P-11 is +0.31 ± 0.05). A metal- 
poor ([Fc/H] = —0.3) isochrone for a = 1.0 from Baraffe et al. 
(1997) is over-plotted as dashed line to investigate the shift due 
to metallicity. Approximately the opposite shift can be expected 
when increasing the metallicity to +0.31, resulting in a mis-match 
between the observational values and the isochrones. 



(1999) Eq. 236 or Tingley & Sackett (2005): 
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Analysis of the transit light curve yields - among other 
parameters - the quantity C/R* that is related to the 
time spent in between the planetary center crossing the 
limb of the star as T& ur — (2^/i?*) -1 . For circular or- 
bits and equatorial transits a/R+ is a function of C/R* 
and the orbital period, and thus can be determined in 
a straightforward way. For eccentric orbits, the relation 
between Q/R* and a/R+, as based on Murray & Dermott 



Thus, in addition to the light curve parameters, a/R+ 
also depends on the orbital eccentricity and argument of 
pericenter, and uncertainties in these parameters propa- 
gate into the error of a/R+. 

In the case of HAT-P-11, the RV data show a sig- 
nificant eccentricity of e = 0.198 ± 0.046 (§ 5). Al- 
though C/i?* and the Vl — b 2 term have small uncer- 
tainties (1.0% and 4.8%, respectively), the uncertainty 
in the (1 + h)/VT 
value of a/R± 

credited to the uncertainties in the orbital parameters, as 
caused by the small RV amplitude and the stellar jitter. 
For comparison, if a/i?* is calculated backwards from the 



e 2 term is higher, yielding a final 
14.6±J;;J. The significant error can be 
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Hipparcos parallax based stellar evolutionary modeling, 
we get a/R-k = 15.58^g2. It is reassuring that this is 
consistent with that derived from the global modeling of 
the data, but the uncertainties are ~ 3 — 4 times smaller. 
This justifies the choice of the Hipparcos-parallax based 
luminosity indicator over using the a/i?* constraint from 
global modeling of the data. This is also confirmed by 
comparing Fig. 5 panels b and c, where the confidence 
ellipsoids from the parallax constraint (panel c) impose 
a tighter constraint on the isochrones. 

3.3. Rotation and Activity of HAT-P-11 

The EPD (prc-TFA) HATNet light curve of HAT-P-11 
shows significant periodic variations with P w 29.2 days 
and a peak-to-peak amplitude of 6.4 mmag in /-band 
(Fig. 7). This period is detected both in the auto- 
correlation function and the Fourier power spectrum of 
the light curve. The period is detected in both seasons 
covered by the HATNet light curve, and the signal re- 
mains in phase across both seasons as well. While this 
period is suspiciously close to the lunar cycle, we note 
that the light curve is relatively unchanged by the EPD 
procedure, i.e. these variations do not appear to correlate 
with any external parameters, including the sky back- 
ground. The TFA procedure does suppress this signal 
when several hundred template stars are used. However, 
this is often the case when applying non-reconstructive 
TFA to other long-period variable stars such as Cepheids 
or Miras. To test whether or not this signal is due to some 
non-astrophysical systematic variation we searched for 
other light curves exhibiting a strong peak in the Fourier 
power spectrum near / = 0.034 [c/d]. We found that 
9 out of the 5000 brightest stars with I < 10.0 showed 
such a peak, however none of these stars are in phase 
with HAT-P-11. Moreover, the distribution of peak fre- 
quencies does not show a pile-up at / = 1/P = 0.034 d _1 
relative to other frequencies, which one might expect to 
see if this variation were a systematic trend. As a final 
test we attempted to recover the signal after applying dif- 
ferent TFA template sets. We tried 100 disjoint template 
sets of 140 light curves. In all cases the P w 29.2 day 
signal was recovered in the auto-correlation function. We 
conclude that the signal is not a systematic variation that 
is present in the light curves of many other stars, and if 
it is not of astrophysical origin then it is due to a phe- 
nomenon that affects the light curve of HAT-P-11 in an 
apparently unique manner. 

A likely interpretation of the variation seen in Fig. 7 is 
that it is due to the rotational modulation of starspots 
on the surface of HAT-P-11. The 6.4 mmag amplitude 
of the variation is consistent with other observations of 
spotted K dwarf stars. Note that there is no signifi- 
cant fi = 0.034 [c/d] periodicity in the RV data, when 
it is modeled as the combination of a Kcplcrian orbit 
plus a sinusoid component (see Fig. 3). The secondary 
peaks in the autocorrelation function in the figure, plus 
the co-phasing of the 29 day variation over two observ- 
ing seasons (left curve of Fig. 7) indicate that individual 
starspots or starspot groups persist for at least several 
rotations. 

If the variation is indeed due to rotational modula- 
tion by starspots, then the rotation period of HAT-P- 
11 is P w 29.2 days. This may be compared with the 
rotation period predicted from the B — V color of the 



star and the Ca II emission index S, using relations of 
Noyes et al. (1984). From Tab. 3, the median value of 
S observed from HIRES is (S) = 0.61. This value and 
B — V = 1.046 ± 0.024 as derived earlier yield a rotation 
period of P ca ic = 24.2 d. The uncertainty in this calcula- 
tion is difficult to quantify. The relations given in Noyes 
et al. (1984) were based on a theoretically motivated fit 
to the rotation period determined for a number of lower 
main sequence stars from the rotational modulation of 
their chromospheric emission. For the 18 K stars in their 
sample with measured rotation periods P b s , the rms dif- 
ference between P b s and the calculated value P ca ic is 
3.5 d. However, Noyes et al. (1984) pointed out that 
the empirical relation for P ca i c as a function of S and 
B — V becomes rather unreliable for B — V > 1.0, be- 
cause of a paucity of stars with observed rotation periods 
in this color range. Hence the true uncertainty is doubt- 
less somewhat larger than 3.5 d. We conclude that a 
rotation period of 29.2 d is consistent with expectations 
from the star's B — V color and (S) value. 

A 29 day rotation period of HAT-P-11, coupled with its 
radius, implies an equatorial velocity v eq = 1.3 km s . 
Let us assume that the stellar rotational axis is inclined 
at i <~ 90° to the line of sight. This is supported by the 
geometry of our Solar system and the generally similar 
geometries of the stellar systems with transiting planets 
whose projected rotational axis inclinations have been 
measured through the Rossiter-McLaughlin effect. We 
would then expect vsini ~ 1.3kms _1 , consistent with 
the value of wsinz reported in Tab. 4. 

Fig. 7 shows the time history of the Ca II S index. 
Fig. 8 shows the very prominent emission cores of the Ca 
II H and K lines observed at three levels of activity during 
the time of the HIRES observations. A long-term varia- 
tion of S is apparent, with timescale close to the length 
of the current data set. As discussed later, this may be 
due to long-term variations of stellar activity analogous 
to the solar activity cycle. Phasing of the S data at a 
period of 29.2 days shows no evidence for periodic behav- 
ior. If the star is really rotating at a period of 29 days, 
this suggests that the chromospheric emission in the H 
and K lines is spread nearly uniformly in longitude over 
the star. 

The color B — V = 1.046 of HAT-P-11 and its me- 
dian level of chromospheric emission (S) = 0.61 imply 
a chromospheric emission ratio R' HK given by log R' HK 
= -4.584 (Noyes et al. 1984). From this we may crudely 
estimate an age for the star using the inverse square root 
relation between activity and age as originally posited by 
Skumanich (1972). Using a fit of this relation to the Sun, 
the Hyades, the Ursa Major group, and 412 individual 
lower main sequence stars as derived by Soderblom, Dun- 
can, & Johnson (1991), we determine a "chromospheric 
age" T cr ^1.25 Gyr. The uncertainty in this estimate is 
difficult to quantify, especially because most of the stars 
in the sample had B — V < 1.0, but it does suggest 
that based on its chromospheric emission level the star 
is likely to be at the low end of the age range given in 
the isochrone fit discussed in § 3.1. 

4. EXCLUDING BLEND SCENARIOS 

4.1. Spectral line-bisector analysis 

As always in determining whether the signature from 
combined photometric and radial velocity variations in 
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Fig. 7.— (Left) The pre-TFA, post-EPD HATnet /-band light curve of HAT-P-11 phased at a period of P = 29.2 days. Only out-of-transit 
points are included in this plot. The apparent 6.4mmag peak-to-peak variation may be due to the rotational modulation of starspots on 
the surface of HAT-P-11. (Upper Right) The auto-correlation function of the light curve at left. Note the first peak at a time-lag of 29 days. 
(Middle Right) The Lomb-Scargle periodogram of the light curve at left, Note the peak at a frequency of 0.034 days -1 . (Bottom) The 
S-indcx as a function of time. Note the long-term variations in the stellar activity. 

sibility that the measured radial velocity variations are 
caused by distortions in the spectral line profiles due to a 
nearby unresolved faint eclipsing binary, whose relatively 
large radial velocity variations mix with the non- varying 
spectrum of the primary and give rise to an apparent ra- 
dial velocity signal with the observed small amplitude. In 
this case it has been shown (e.g. Queloz et al. 2001) that 
the bisector span (BS) of spectral lines in the blended 
spectrum varies in phase with the radial velocity signal 
itself, with similar amplitude. We have carried out an 
analysis of the bisector span based on the Keck spec- 
tra as described in earlier HATNet detection papers (see 
Bakos et al. 2007a), and do not see any statistically sig- 
nificant correlation between the line bisector spans and 
the measured radial velocities (confer Fig. 4 and Tab. 3), 
thus providing no support for the hypothesis that the 
signal is caused by a blended eclipsing stellar system. 

However, because of the small amplitude of the RV 
signal (and hence small amplitude of expected bisector 
span variations if the apparent RV signal is indeed due 
to a blend), coupled with the large jitter as seen in the 
RV residuals in Fig. 4, the above test cannot completely 
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Fig. 8. — The Calcium K (left) and H (right) line profile in se- 
lected HIRES observations of HAT-P-11. Both panels show three 
spectra overlaid; data taken at high, median and low activity, as 
characterized by the S index. The spectra are matched to a com- 
mon flux/ wavelength scale using points outside the H and K line 
cores. The vertical axes on the plots are arbitrary and proportional 
to counts. 

a star is due to a transiting planet, it is necessary to 
exclude the possibility that the entire combined infor- 
mation is due to a set of circumstances that falsely give 
rise to the characteristic signature of a transiting planet. 
Following Torres et al. (2007), we first explored the pos- 
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rule out contamination by an unresolved binary system. 
Thus we sought other, independent ways of excluding 
blend scenarios. In the following subsections we attempt 
to model the system as a hierarchical triple, where the 
smallest component is either a star or a Jupiter-sized 
planet, and also as a background eclipsing binary blended 
with the light of a foreground K dwarf. We show that 
none of these models are consistent with all of the avail- 
able observations. 

4.2. Detailed blend modeling of a hierarchical triple 

We consider the possibility that HAT- P- 11 is a hier- 
archical triple system, having a deeper intrinsic eclipse 
of two bodies diluted by the light of the bright K dwarf. 
Note that the bright K dwarf with well known proper- 
ties (parallax, proper motion, etc.) in this putative triple 
system is referred to as HAT-P-11 in the following dis- 
cussion. To rule out the hierarchical triple scenario, we 
attempt to fit a blend model to the observations. There 
are two scenarios for the eclipsing system in the triple 
system. The first one is a stellar eclipsing binary. The 
second case is a transiting hot Jupiter orbiting a low- 
mass star, with the few percent deep transit diluted by 
HAT-P-11 to become only 4.2mmag. We note that the 
shallower the transit event, the larger the variety of con- 
figurations that can match the observations. HAT-P-11 
represents a transition towards millimag transit events, 
the blend analysis of which will certainly be challenging. 

In each case (stellar binary, vs. transiting hot Jupiter 
blend) we fit the follow-up z, r and ip-band light curves 
together with the HATNet /-band light curve following 
a procedure similar to that described by Torres et al. 
(2005). We assume that the bright star (HAT-P-11), 
which is not eclipsed, has the mass, metallicity, age and 
distance determined in § 3. We also assume the compo- 
nents of the eclipsing system have the same metallicity, 
age and distance as the bright star, i.e. they form a hier- 
archical triple. We fix the ephemeris, orbital eccentricity 
and argument of periastron to the values determined in 
§ 5.4 and vary the masses of the two components and the 
inclination of the orbit. The downhill simplex algorithm 
is used to optimize the free parameters. For the case 
where one of the components is a planet, we vary the 
radius of the planet rather than its mass. For the eclips- 
ing binary case the magnitudes and radii of the stars are 
taken from the Padova isochrones (Girardi et al. 2000) 
applying the radius correction described by Torres et al. 
(2005). We use these isochrones to allow for stars smaller 
than 0.4 M©. For the planet case we use the Yonsei-Yale 
isochrones to be consistent with the single-star planet 
modeling. 

Fig. 9 shows the best fit model for the eclipsing bi- 
nary case and two illustrative models for the blended 
hot-Jupiter case over-plotted on the z-band light curve. 
We are unable to fit the light curve using a combina- 
tion of three, physically associated stars, but are able 
fit it as a planet transiting one component of a binary 
star system, so long as the star hosting the planet has 
M > 0.72 M Q . Interestingly, 0.72 M© emerges as a crit- 
ical mass (and corresponding critical radius, luminosity) 
for the eclipsed star, and splits the parameter space into 
two disjoint domains. As the mass of the host star is 
reduced below 0.72 M© (and its radius is decreased), the 
only way to maintain the depth and duration of the al- 



ready central (b = 0, i.e. longest possible) transit is to 
increase the radius of the planet. This, however results 
in a longer ingress and egress time than is allowed by the 
observations. Above 0.72 M© the radius of the planet 
needed to fit the transit depth changes slowly, while rais- 
ing the impact parameter of the planet can compensate 
to fit the transit duration (see Fig. 9). While a model 
of this form fits the light curve, the F-band light ratio 
of the planet-hosting star to the brighter 0.81 M© star 
would be > 0.46. The spectrum of a second star this 
bright should have been easily detected in the Keck and 
DS observations, and therefore we can rule out both of 
these blend scenarios. 
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Fig. 9.— (Top) The z-band follow-up light curve for HAT-P-11 
with model light curves for three blend models over-plotted. The 
dark points show the binned light curve , while the light points 
show the unbinned light curve. The solid line shows a model where 
a planet orbits one component of a binary star, with the planet- 
hosting star having a mass of 0.8 Mq. The dot-dashed line shows 
a similar scenario but where the planet-hosting star has a mass 
of 0.4 Mq . The dashed line shows the best fit hierarchical triple 
model consisting of eclipsing binary star system bound to a distant 
third star. We can rule out the hierarchical triple scenario and the 
blend models where a planet orbits a binary star component that is 
smaller than 0.72 Mq, blend scenarios with a larger planet-hosting 
star can fit the light curve but a second star this massive would 
have been detected in the Keck and DS spectra. (Bottom) Here 
we show x i the planet radius R p and the impact parameter b as 
a function of the planet-hosting star mass for the blend scenario 
of a planet orbiting one component of a binary system. \ 2 ' s 
linearly scaled between 95788 and 100457 while R v is linearly scaled 
between 0.583-Rj and 2.17-Rj. For stellar masses below 0.72 Mq 
the planet radius needed to fit the transit depth and duration yields 
an ingress and egress time that is too long. Above 0.72 Mq the 
impact parameter increases to fit the transit duration while the 
planet radius needed to fit the transit depth changes slowly. 



4.3. Contamination from a background eclipsing binary 

An alternate model of an astrophysical false positive 
that would make the contaminating star much fainter rel- 
ative to the main object is one involving a background 
eclipsing binary (chance alignment). We modeled this 
case using the same methodology described above, fit- 
ting the combined z-band follow-up light curves as the 
sum of the light from three stars. The properties of the 
main star were held fixed as before, and those of the 
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Fig. 10. — Sample blend model fitted to our combined 2-band 
measurements (arbitrary scale on the vertical axis). The solid line 
represents the light curve resulting from a background eclipsing 
binary system whose flux is diluted by the main star, which con- 
tributes 95% of the 2-band light in this case. The eclipsing binary 
is composed of 0.64 Mq and 0.13 Mq stars in an edge-on config- 
uration, placed 60 pc behind the target. The predicted relative 
brightness of the primary in this scenario is only 3% of the light of 
the main star in the optical. 

eclipsing binary components were constrained to lie on 
the same isochrone, which for simplicity was taken to 
be the same as the main star. Extensive tests indicated 
that the detailed shape of the resulting synthetic light 
curve is fairly insensitive to the distance at which we 
place the eclipsing binary behind the target, and thus 
acceptably good fits to our z-band photometry can be 
achieved for a wide range of distances, given the mea- 
surement precision. Fig. 10 shows an example of such a 
fit, in which the eclipsing binary is placed about 60 pc 
behind the main star and is composed of a 0.64 Mq star 
(spectral type K7V-M0V) orbited by a 0.13 M stellar 
companion. This model is not only consistent with the 
photometry, but it predicts an optical brightness for the 
primary of the eclipsing binary of only 3% relative to the 
main star. Detecting such a faint set of spectral lines in 
our Keck spectra would be challenging. Furthermore, if 
we place the eclipsing pair 110 pc behind the main star, 
the light curve fit is still about the same, but the rela- 
tive brightness decreases by a factor of two, making the 
eclipsing binary spectroscopically undetectable. Spectral 
line bisector variations predicted by this second model 
(see, e.g., Torres et al. 2005) would be at the level of the 
scatter in our actual measurements, and thus could not 
be entirely ruled out either. 

While this blend scenario appears to satisfy all observa- 
tions, it implies the presence of an eclipsing binary ^4.5 
magnitudes fainter in the optical very near the present 
location of our target. It is fortunate that HAT-P-11 has 
a large proper motion(0'/264 yr _1 ; Perryman et al. 1997). 
Using Palomar Observatory Sky Survey plates from 1951 
(POSS-I, red and blue plates), we can view the sky at the 
current position of HAT-P-11 unobstructed, since the tar- 
get was 15" away. Fig. 11 shows a stamp from the POSS- 
I, POSS-II plates, and also a current observation with 
the FLWO 1.2 m telescope. The proper motion of HAT- 
P-11 is apparent. A number of faint stars are marked on 
the figure. Note that the faintest possible star that can 
cause the observed 4.2 mmag variation in the combined 
light curve would be 6 magnitudes fainter than HAT-P- 
11, and in this extreme situation the blended faint star 
would need to completely disappear during its transit. 



There is no faint star down to ~ 19 mag within ~ 5" 
of the current position of HAT-P-11. The closest star 
is marked "2" in Fig. 11 is 2MASS 19505049+4805017, 
with z w 14.4, or about z w 5.6 magnitudes fainter 
than HAT-P-11. This is well resolved in some of our 
follow-up photometry observations, such as on the 2007 
September 2 night, and its brightness was constant, with 
r.m.s. much smaller than the required amplitude of its 
variation should be (<~ 1.3 mag) to cause the observed 
4.2 millimag combined dip. 




Fig. 11.— Images of a 2' X 1.7' field containing HAT-P-11 from 
the POSS-I Red survey (Left), POSS-II Red survey (Center) and 
a FLWO 1.2 m z-band image (Right, see § 2.3). The dates of the 
exposures are 10 July 1951, 25 August 1989, and 2 September 
2007 respectively. The cross marks the position of HAT-P-11 in 
the POSS-I image (labeled as star 1), while the diamond marks 
the position of HAT-P-11 in the POSS-II image. Between 1951 and 
2007 HTR155-001 moved 13.8" in the NNE direction. Stars labeled 
1 through 4 on the POSS-I image have USNO A-2.0 R magnitudes 
of 10.8, 15.1, 17.2 and 18.9 respectively (USNO-A2.0; Monet 1998). 
From the POSS-I image we can rule out stars brighter than R ~ 19 
at the current position of HAT-P-11. 

In summary, i) we have not seen bisector variations 
correlated with the radial velocity, ii) we investigated 
unrcsolvable hierarchical triple systems (blended eclips- 
ing binary or transiting hot Jupiter), and found all to 
be incompatible with the photometric and spectroscopic 
data together, and iii) we were able to exclude chance 
alignment of a background eclipsing binary based on the 
high proper motion of the star. We therefore conclude 
that the transit signal, and the synchronized RV signal, 
indeed are both due to a sub-stellar companion transiting 
the K4 dwarf HAT-P-11. 

5. GLOBAL MODELING OF THE DATA 

In order to perform an optimal analysis of the data at 
hand, we assembled a global model of a planetary transit 
scenario that describes all the data components: i) the 
HATNet discovery light curve, ii) the photometric follow- 
up light curves of 10 independent events, and iii) the ra- 
dial velocity data from Keck (§ 5.1). The physical model 
was then extended by a model that describes systematic 
variations in the data (§ 5.2). We then performed a joint 
fit of the combined physical plus systematic model by 
using all these data simultaneously to derive the phys- 
ical parameters, and also to correct for the systematic 
variations (§ 5.3). The results are discussed in § 5.4 and 
§ 5.5. 

5.1. Physical model for the light curve and radial 
velocity 

Our model for the follow-up light curves used analytic 
formulae based on Mandel & Agol (2002) for the eclipse 
of a star by a planet, where the stellar flux is described 
by quadratic limb-darkening. We did not attempt to 
fit for the limb darkening coefficients, as they require 
very high quality data to resolve their degeneracy with 
other parameters. Instead, the limb darkening coeffi- 
cients were derived from the second iteration SME results 
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(§ 3), using the tables provided by Claret (2004) for z, r 
(FLWO 1.2m) and I c bands (Konkoly 0.6m Schmidt). 
The transit shape was parametrized by the normalized 
planetary radius p = R p /R+, the square of the impact pa- 
rameter b 2 , and the reciprocal of the half duration of the 
transit 16 (,/R*. We chose these parameters because of 
their simple geometric meanings and the fact that these 
show negligible correlations (see Carter et al. 2008; Bakos 
et al. 2007b; Pal 2008b). 

Our model for the HATNet data was a simplified ver- 
sion of the Mandel & Agol (2002) analytic functions, be- 
cause i) the number of in-transit data points in the HAT- 
Net light curve is significantly less than the same num- 
ber for the follow-up light curves, and ii) the individual 
errors on the points is much worse for HATNet than for 
the follow-up. Thus, while the HATNet data are efficient 
in constraining the cphemcris, they are not comparable 
to the follow-up data in determining other transit pa- 
rameters that depend on the exact shape of the transit 
(such as depth, or duration of ingress). For these rea- 
sons, we employed a simple transit model that neglects 
limb darkening. This model function, which we label 
as "P1P3", originates from an expansion of the exact 
Mandel & Agol (2002) model function by Legendre poly- 
nomials, resulting in a functional form (omitting scaling 
factors) of the ingress (and symmetrically, the egress) of 
x — > (21x — 5x 3 )/16, where the independent variable x is 
scaled to be 1 at first (third) contact, and —1 at second 
(fourth) contact. For reference, the trapeze model has a 
simple linear function describing the ingress: x — ► x. The 
"P1P3" model approximates the Mandel & Agol (2002) 
function to better than 1%, and our experience shows 
that such an approximation yields the same timing preci- 
sion as more sophisticated model functions. In addition, 
using the simplified model yields much faster computa- 
tion. A similar transit model approximation using hyper- 
bolic tangent functions was described earlier by Protopa- 
pas, Jimenez, & Alcock (2005). The depth of the HAT- 
Net transits was also adjusted in the fit (independent of 
the follow-up data) and we found that the depth derived 
from the HATNet light curve is comparable (§ 5.4) to 
the depth of the (EPD and TFA reconstructed) follow- 
up transit measurements, namely the ratio of the two 
was Binst = 0.95 ± 0.08. In general, this "instrumen- 
tal blend" parameter has to be adjusted independently 
since the possible contamination by nearby stars in the 
wide-field images may yield shallower transits. 

The adopted model for the radial velocity (RV) varia- 
tions was described earlier in § 2.4. The RV curve was 
parametrized by an eccentric Keplerian orbit with semi- 
amplitude K, RV zero-point 7, Lagrangian orbital ele- 
ments (k, h) = e x (coscj, sinu>), plus a linear trend G\. 

We assumed that there is a strict periodicity in the 
individual transit times. In practice, we introduced the 
first transit center as Ta — T" Ci _23i and the last tran- 
sit center as Tb = T Cj+ gi, covering all of our measure- 
ments with the HATNet telescopes, the FLWO 1.2 m 
telescope and the Schmidt telescope. The transit cen- 
ter times for the intermediate transits were interpolated 
using these two epochs and the N tr transit number of 

16 The duration is defined as the time between the contact cen- 
ters, i.e. when the center of the planetary disk crosses the limb of 
the star during ingress and egress. 



the actual event under the assumption of no transit tim- 
ing variations. The model for the RV data contained the 
ephemeris information through the T Ci _23i and T c<+ gi 
variables; for instance the period assumed during the fit 
was P = (T C)+9 i - T c _23i)/(91 + 231). The other cou- 
pling between the transit photometry and the RV data 
is via the k and h Lagrangian orbital elements that de- 
termine the relation between the photometric and or- 
bital ephemeris, plus have a minor effect on the tran- 
sit shape. Altogether, the 12 parameters describing the 
physical model were T Ci _ 2 3i, T c>+9 i, R p /R+, b 2 , C/-R*, 
K , 7, k — ecoscj, h — esinw and G\. As mentioned 
above, the parameters were extended with the instru- 
mental blend factor B- msi , and the HATNet out-of-transit 
magnitude, M ,HATNot- 

5.2. Models for systematic variations 

A "joint" fit of the physical model to the global data- 
set has been performed routinely for recent discover- 
ies (Bakos et al. 2008; Latham et al. 2008). Here we 
extended our physical model that describes the tran- 
sit events and the radial velocities with an instrumen- 
tal model that describes the systematic variations of the 
data. Since the HATNet photometry has been already 
EPD and TFA corrected (§ 2.1), and'the Keck RVs have 
been investigated for trends (§ 2.4), but deemed to have 
insufficient number to decipher any significant system- 
atic variation (50 points compared to 0(10 ) points for 
photometry), we only modeled systematic variations of 
the follow-up photometry (§ 2.3). As mentioned earlier, 
the effect of systematics on the shallow transit in the 
photometric follow-up data is enhanced, and can impact 
the basic planetary parameters, such as the planetary 
radius. We experimented with several models, since the 
treatment of systematics may improve precision, but not 
necessarily the accuracy of the results. For reference, we 
call the simplest model, with no EPD and no TFA, as 
model "E0T0" (for technical details, see § A). 

Model "E0TRG" was independent (i.e. not simultane- 
ous) reconstructive global TFA of the FLWO 1.2 m z- 
band data (8 nights, see Tab. 1), assuming a trapeze 
shaped mo (p, t) model function in Eq. A2 with three 
parameters (total duration, ingress duration, depth). 
The period was iteratively refined using both the 
FLWO 1.2 m and the HATNet data, but reconstructive 
TFA was only performed on the FLWO photometry. We 
used the magnitude values transformed to the selected 
reference frame for each night (see § 2.3), and the zero- 
point shift for each night was determined separately. We 
increased the number of template stars till the S/N of 
the trapeze reached its maximum (S/N = 39), requiring 
246 stars with z < 14. Note that this was a global TFA, 
in the sense that there was one Cfc TFA coefficient per 
star in Eq. A3 for all nights combined (see § A. 2). The 
reconstructed light curve was then fed into the joint fit. 

Model "ELT0" had EPD parameters in the simulta- 
neous fit, together with the parameters of the physical 
model, but no TFA-based detrending, i.e. all the system- 
atic variations in the follow-up photometry were mod- 
eled as due to external parameters. We chose 5 such pa- 
rameters, namely the hour angle (characterizing a mono- 
tonic trend that linearly changes over time), the square 
of the hour angle, the stellar profile sharpness parameter, 
S = (2.35/FWHM) 2 , the X and Y pixel position on the 
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chip (each with two coefficients). The exact functional 
form of the above parameters contained 8 coefficients, 
including the auxiliary out-of-transit magnitude of the 
individual events. The EPD parameters were indepen- 
dent for all 10 nights, implying 80 additional coefficients 
in the global fit. 

Model "EOTL" had local TFA parameters in the instru- 
mental model for each night, and no EPD was used. We 
selected 23 TFA template stars that were present in all 
observations, representing faint and bright stars spread 
across the chip. The total number of coefficients in 
the simultaneous fit was thus 230. Model "EOTG" had 
global TFA parameters for all nights, using 211 stars that 
were present on all frames. Note that without simulta- 
neous fit using an underlying physical model (nij(ti)) in 
Eq. A2) that takes into account the different shape of 
the light curves acquired through different photometric 
bands, such a global TFA fit would be not possible. 

We tested variants of the above methods. Model 
"ELTL" was a simultaneous fit of EPD and TFA param- 
eters in local EPD-mode and local TFA-mode with alto- 
gether 80 + 230 parameters in addition to the 12 physical 
parameters. Model "ELTG" was a simultaneous fit of lo- 
cal EPD and global TFA parameters with an additional 
80 + 211 instrumental parameters. The parameter sets 
were extended in all cases with the out-of-transit mag- 
nitudes for each follow-up light curve and the HATNet 
photometry. The number of fitted parameters was much 
smaller than the number of data points (~ 5000). 

5.3. Performing the joint fit 

It is computationally rather intensive to perform a fit 
with ~ 300 parameters and to determine their respective 
error distributions. 

Since the majority of the fitted parameters appear as 
linear terms in the final form of the model functions, we 
minimized \ 2 m the parameter space by using a hybrid 
algorithm, namely by combining the downhill simplex 
method (a.k.a. AMOEBA, see Press et al. 1992) with 
the classical linear least squares algorithm. The simplex 
itself was propagating in the "nonlinear" hyperplane of 
the parameter space, while in each point of the hyper- 
plane the value of \ 2 was computed using classical linear 
least squares (CLLS) minimization. This method yielded 
a very good convergence even with the large number of 
free parameters. 

Uncertainties of the parameters were derived using the 
Markov Chain Monte-Carlo method (MCMC, see Ford 
2006), by starting two types of chains from the best fit 
value. The first one was the classical MCMC chain, let- 
ting all parameters vary. The other type of chain was 
generated by only allowing the variation in a hyperplane 
of the parameter space consisting of the nonlinear pa- 
rameters and of some additional linear parameters where 
error-calculation was desired (examples are K, 7 or G\). 
In this "Hyperplane-CLLS" chain the remaining linear 
parameters were determined at each step via CLLS. The 
H-CLLS chain has the advantage that transition proba- 
bilities are much higher in the separated hyperplane, and 
thus the convergence and computing time requirements 
are also smaller by an order of magnitude. Furthermore, 
the error distribution of parameters for the hyperplane 
variables was identical to those determined by classical 
MCMC. The distribution for the parameters outside the 



hyperplane was distorted, but these were considered as 
auxiliary, where such distribution is irrelevant (e.g. the 
EPD and TFA coefficients). Therefore the final distri- 
bution of the physical parameters was derived by using 
H-CLLS chains. 

For both types of the chains the a priori distributions 
of the parameters were chosen from a generic Gaussian 
distribution, with eigenvalues and eigenvectors derived 
from the Fisher covariance matrix for the best fit value. 
All functions that appear in the physical and instrumen- 
tal model are analytical, and their partial derivatives are 
known. The parametric derivatives of the transit light 
curve model of Mandel & Agol (2002) are found in Pal 
(2008b), while the parametric derivatives of the eccen- 
tric radial velocity curves can be obtained by involving 
the implicit function theorem, and can be expressed as 
functions of the eccentric anomaly. Thanks to these an- 
alytic properties, the derivation of the Fisher matrix was 
straightforward. 

5.4. Results of the fit and the effect of systematics 

We performed a joint fit on the HATNet light curve, 
the follow-up photometry, and the Keck radial velocities 
(the "data"), using the physical model described in § 5.1, 
as extended by the various models for systematics (§ 5.2), 
and using the algorithms described in § 5.3. 

We carried out the fit considering each model describ- 
ing the systematics in the follow-up photometry sepa- 
rately in order to investigate their effect on the final re- 
sults. As mentioned earlier, this is important due to the 
very shallow transit signal and the large relative ampli- 
tude of systematic variations. We used the m/ fitted 
magnitude values (§ 2.1) for the photometry follow-up 
light curves (without prior EPD or TFA) . Since only the 
follow-up photometric data-component was described by 
the different systematic models (see § 5.2), we focused on 
those light curve parameters that are affected by the pho- 
tometric follow-up, namely the transit center time of the 
last event T c + gi, the radius ratio of the planet to star, 
the (square of the) impact parameter and the transit du- 
ration parameter £/-R*- The other adjusted parameters 
were not relevant in this comparison since they are not 
affected by the follow-up light curves. Note, that with 
the exception of the E0TRG model, all EPD and TFA pa- 
rameters were fit simultaneously with the physical model 
parameters. 

The best fit values and the respective uncertainties 
are summarized in Tab. 5. It is clear that EPD and 
TFA both significantly reduce the uncertainties in all 
the parameters. Note that the implementation of E0TRG 
was somewhat sub-optimal. Because the reconstructive 
trapeze shape was not bandpass-dependent (as opposed 
to the full Mandel & Agol (2002) model used in the other 
cases), only a fraction of the follow-up was used, namely 
all photometry in z-band. When using EPD and TFA 
together, the global TFA (ELTG) performs better than 
local TFA (ELTL) for each night. The improvements by 
EPD (ELTO) and TFA (EOTG) separately are roughly the 
same. The best model, based on the formal error-bars, is 
the EPD and global TFA together (ELTG); here the un- 
biased error in the light curve parameters decreased by 
a factor of 1.5. 

The final value for the parameters also changed by 3-a 
in the sense that both i? p /-R* and b 2 decreased. This 
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TABLE 5 

LIGHT CURVE PARAMETERS INVOLVING DIFFERENT KINDS OF MODELS FOR 
SYSTEMATIC VARIATIONS. 



Method T Cj+9 i - 2454790.0 R P /R* b 2 C/R* 

E0T0 a 1.62781 ± 0.00055 0.06035 ± 0.00140 0.165 ± 0.120 22.05 ± 0.24 

E0TL b 1.62793 ± 0.00047 0.05831 ± 0.00094 0.114 ± 0.087 21.82 ±0.18 

E0TG c 1.62772 ± 0.00046 0.05730 ± 0.00096 0.112 ± 0.091 21.99 ±0.19 

E0TRG d 1.62758 ± 0.00133 0.05690 ± 0.00150 0.223 ± 0.137 22.10 ±0.27 

ELT0 e 1.62834 ± 0.00044 0.05825 ± 0.00094 0.109 ± 0.093 22.25 ±0.17 

ELTL f 1.62826 ± 0.00040 0.05675 ± 0.00097 0.113 ± 0.093 22.17 ±0.16 

ELTG g 1.62832 ± 0.00039 0.05758 ± 0.00091 0.120 ± 0.087 22.15 ±0.15 



a No EPD and no TFA performed on the photometric follow-up data. 

b No EPD, and per-night "local" TFA performed with 23 TFA template stars and 
simultaneous fitting with the transit model. 

c No EPD, and global simultaneous TFA performed with 211 TFA templates. 

d No EPD. Global reconstructive trapeze TFA performed with 246 TFA templates. 
Only the 8 nights of FLWO z-band observations were used. 

C EPD performed using 5 free parameters per night. No TFA. 

f EPD as for ELT0. TFA as for E0TL. 

SEPD as for ELT0. TFA as for E0TG. 



means that the treatment of the systematics recovers 
a sharper transit with smaller impact parameter and 
smaller planetary radius. If the stellar parameters were 
determined using the a/R+ luminosity constraint (§ 3.2), 
then ignoring the correction for systematics would lead 
to higher impact parameter, larger R p /R+, and smaller 
a/R+, corresponding to smaller density (Eq. 1). One may 
expect that the effect of systematics is slowly diminished 
by accumulating more data taken with different instru- 
ments. This is consistent with our experience, that is the 
initial follow-up light curves without proper treatment 
of systematics indicated larger impact parameters, and a 
lower density dwarf star that was inconsistent with the 
Hipparcos-based stellar parameters (§ 3.1). By accumu- 
lating more data this inconsistency diminished (§ 3.2). 
Because both EPD and TFA point toward this direction, 
and also deliver the smallest formal uncertainty, we ex- 
pect that the most accurate model will be the ELTG in 
Tab. 5. This conjecture will be decided with accurate 
ground-based follow-up photometry (such as presented 
recently in Johnson et al. 2008a), and ultimately when 
the Kepler space mission returns the light curve for HAT- 
P-ll. This will have important implications on optimal 
trend removal for other transiting systems, where such a 
space-born "reference" will not be available. 

The results for the simultaneous fit using the ELTG sys- 
tematic model for the follow-up light curves were the 
following: T c ._23i = 2453217.75466 ± 0.00187 (BJD), 
T Cj+9 i = 2454605.89132 ± 0.00032 (BJD), K = 11.6 ± 
1.2ms -1 , k = ecosw = 0.201 ± 0.049, h = esinw = 
0.051 ± 0.092, R p /R+ = 0.0576 ± 0.0009, b 2 = 0.120 ± 
0.087, C/#* = 22.15 ± 0.15day-\ and 7 = -0.4 ± 
0.8ms' 1 , B ins tr = 0.95 ± 0.08, Gi = 0.0297 ± 
0.0050 ms-7day, M 0jH ATNct = 8.35892 ± 0.00003 (in- 
strumental HATNet out-of-transit magnitude). The 
combined and phase-binned light curve of all z-band 
FLWO 1.2 m observations is shown on Fig. 12. The plan- 
etary parameters and their uncertainties can be derived 
by the direct combination of the a posteriori distribu- 
tions of the light curve, radial velocity and stellar param- 
eters (see also Pal et al. 2008a). We found that the mass 
of the planet is M p = 0.081 ± 0.009 M 3 = 25.8 ± 2.9 M e , 
the radius is R p = 0.422 ± 0.014 R 3 = 4.73 ±0.16^0 and 



its density is p p — 1.33±0.20gcm 3 . The final planetary 
parameters are summarized at the bottom of Table 6. 
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Fig. 12. — Combined and binned light curve of HAT-P-11, in- 
volving all the z band photometry, obtained by KeplerCam. The 
light curve is superimposed with the best fit model (which is the 
result the joint fit described in § 5). 



5.5. Constraints on orbital eccentricity 

An interesting aspect of the HAT-P-11 system is that 
the light curve analysis, the Hipparcos parallax and the 
theoretical isochrones together provide an extra con- 
straint on the orbital eccentricity. In our analysis the or- 
bital eccentricity and argument of pericenter are charac- 
terized by the Lagrangian orbital elements k — e cos u> — 
0.201 ±0.049 and h = esinw = 0.051 ±0.092. If the anal- 
ysis relied purely on the radial velocity data, then the 
errors in k and h would be similar, the error of h being 
somewhat smaller (see § 2.4). However, when the data 
are complemented with photometry of transit events (as 
in our case), the period and phase of the orbit become 
tightly constrained. As a result, the k and h orbital 
elements will have different uncertainties since the mean 
longitude at the transit is constrained more by cos ui than 
sin u (this is due to the same underlying reason that the 
phase lag between the secondary and primary transits is 
also proportional to e cos u> and does not strongly depend 
on esinw). Indeed, the a posteriori distribution of h and 
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TABLE 6 
Orbital and planetary parameters 



Parameter 



Value 



Light curve parameters 

P (days) 4.8878162 ± 0.0000071 

T c (BJD) 2454605.89132 ± 0.00032 

T M (days) a 0.0957 ± 0.0012 

T12 = T 34 (days) a 0.0051 ± 0.0013 

'' '<'■ 15.58tg;g 

C/R* 22.15 ±0.15 

Rp/R+ 0.0576 ± 0.0009 

b = acos i/R+ 0.347lo 

i (deg) 88.5±o'.6 

Spectroscopic parameters 

K (ms- 1 ) 11.6 ± 1.2 

7 (kms -1 ) -0.4 ±0.8 

Gi (ms _1 /day) 0.0297 ± 0.0050 

fc RV b 0.201 ± 0.049 

/i RV b 0.051 ± 0.092 

k c c 0.190 ±0.046 

h c c -0.016 ± 0.056 

e 0.198 ±0.046 

uj 355.2 ± 17.3° 

Secondary eclipse parameters' 1 

T s (BJD) 2454608.96 ± 0.15 

T a 14 0.1006 ± 0.0130 

T M2 0.0054 ± 0.0013 

Planetary parameters 

M p (Mj) 0.081 ± 0.009 

Rp (Rj) 0.422 ±0.014 

C(Mp,R p ) d 0.025 

p p (gem" 3 ) 1.33 ±0.20 

" (AU) 0.05301°;°°°! 

logg p (cgs) 3.05 ±0.06 

T eq (K) 878 ± 15 

e 0.025 ±0.003 

F per (ergs" 1 cm" 2 ) ° 2.04 • 10 s ± 2.78 • 10 7 

F ap (ergs- 1 cm- 2 ) ° 9.11 • 10 7 ± 9.53 • 10 6 

(F) (ergs- 1 cm" 2 ) c 1.34 • 10 s ± 9.38 • 10 6 

a Ti4: total transit duration, time between first to last contact; 
T12 = T34: ingress/egress time, time between first and second, 
or third and fourth contact. 

b Lagrangian orbital elements, based purely on the RV data. 

c Refined value of h and k, derived from RV analysis and the 
constraint given by the Ckh = (Vl — e 2 / (1 + h) value, resulting 
from the light curve modeling and stellar evolution analysis. 

d Correlation coefficient between the planetary mass M p and 
radius Rp. 

e Occultation parameters and incoming flux per unit surface 
area in periastron, apastron, and average for the orbit were 
calculated using the refined values kc and he- 

k as yielded by the joint analysis (§ 5.4) is asymmetric, 
and the uncertainty in k is half of the uncertainty in h 
(Fig. 13, large oval around the filled diamond). 
By re-arranging Eq. 2, we can define Ckh as: 



Ckh 



Vl-h 2 - k 2 
T+h 



A (A. 



(3) 



where n = lixjP is the mean motion (see also Ford, 
Quinn, & Vcras 2008). The period P, impact param- 
eter b, and are well constrained by the photometry 
of multiple transit events. The normalized semi-major 
axis a/R+ is determined from the apparent magnitudes, 
Hipparcos parallax and isochrones, as described in § 3.2. 
These together yield an estimate of Ckh = 0.931 ± 0.063, 
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Fig. 13. — The cloud of points represent the a posteriori distri- 
bution of the orbital elements k = ecosw and h = esinw yielded 
by the global modeling (§ 5.4). The best fit (k, h) pair is marked as 
a filled diamond, encircled with the 1-cr confidence ellipsoid. The 
constraint given by the Ckh = Vl — e 2 /(l + h) quantity defines 
a stripe in the (h, k) plane. The location of the finally accepted 
(fc, h) value and its 1-cr confidence are denoted by the dot below 
the diamond, and the surrounding ellipse, respectively. 



independent of the RV data, and thus provide an addi- 
tional constraint on h and k in the form of a hyperbola 
on the k — h plane. Taking into account the uncertain- 
ties in Ckh, this appears as a stripe in the k — h plane 
(Fig. 13). This constraint can be used to decrease the 
uncertainty of h by a factor of two, so that it becomes 
similar to that of k. 

In practice, we have a Monte Carlo distribution for k, 
h and Ckh- For each k and h value in MC distribution, 
we determine the closest k' ,h! point of the hyperbola 
(there is a unique solution for this). Since the value for 
Ckh constrains only one degree of freedom, we can treat 
(fc, h) and (fc', h') as independent. Therefore, we define 
kc and he as the mean of k, kl and h, h', respectively, 
and consider these as an improved estimate of the orbital 
eccentricity parameters. The result is an MC distribution 
for kc and he, where their mean values and errors can 
be derived: k c = 0.190±0.046 and h c = -0.016±0.056. 
Indeed, the errors in kc and he are now comparable. 

Based on kc and he, the refined orbital eccentricity 
is e = 0.198 ± 0.046 and the argument of pericenter is 
uj = 355. 2± 17.3°. Note that the above method is unique, 
and can be applied for eccentric transiting systems with 
well-known parallax but poorly determined RV orbit. It 
is likely that the Kepler mission will find similar cases. 
A counter-example is HAT-P-2b (Bakos et al. 2007b), 
where the RV amplitude is large compared to the uncer- 
tainties (1 kms -1 ), and the error in the parallax is larger, 
thus h can not be refined. The current analysis points 
towards an even more general way of simultaneously fit- 
ting to all data, including the stellar isochrone search 
with various constraints (e.g. parallax information with 
errorbars). 
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The stellar parameter determination (§3), the blend 
analysis (§ 4) and the global modeling of the data pre- 
sented in § 5 together show the detection of a 25.8 ± 
2.9 M e and 4.73 ± 0.16 i?© planet that orbits a bright 
(V = 9.59) K4 dwarf star on an eccentric (e — 0.198 ± 
0.046) orbit with P — 4.8878 day period, causing a 
11.6 ± 1.2ms _1 RV variation of the host star and a 
4.2mmag transit as it passes in front of the star. The 
planet is somewhat more massive than our own Uranus 
(14.3 M e ) and Neptune (17.1 M e ), and slightly more 
massive than GJ 436b (22.6 Mq). Based on its mass and 
radius, it is justified to use the metonym super-Neptune 
for classification of HAT-P-llb, just like for GJ 436b, 
the only similar object known so far. For the location of 
~Neptune-mass objects, see Fig. 14. 

When compared to models of Fortney et al. (2007), 
HAT-P-llb is much smaller in radius than similar mass 
planets with 50% rock/ice core and 50% H/He envelope, 
and HAT-P-llb is much bigger than pure rock/ice plan- 
ets without a gas envelope. The mismatch with such 
pure- rock/ice "super-Earths" is confirmed by models of 
Valencia, Sasselov, & O'Connell (2007). Overlaid in 
Fig. 14 are super-Earth planet models with Earth-like 
composition (lower dotted line) and 50% by mass H 2 
content and Earth-like Si/Fe ratio (upper dotted line). 
Both models fall clearly below the observational values 
for HAT-P-llb. As the expected limit for gravitational 
capture of H/He is ~ 10M© (Rafikov 2006), it is hardly 
surprising that HAT-P-llb has some H/He, even if evap- 
oration played an important role during its evolution. 
HAT-P-llb receives only modest irradiation when com- 
pared to hot Jupiters. The time-averaged flux over its 
eccentric orbit is 1.34 • 10 8 ± 9.38 • 10 6 ergs -1 cm~ 2 , and 
the equivalent semi-major axis where similar irradiation 
would be received from a Sun-like star is a equ i V ~ 0.1 AU. 
Using the models from Fortney et al. (2007), the ra- 
dius of HAT-P-llb is indeed consistent with a 90% 
heavy element planet with an irradiation corresponding 
to a equiv — 0.1 AU. In these models, the heavy elements 
are located in the core, and the envelope is metal-free. 

We also compared the observational values to the thor- 
ough theoretical work of Baraffe et al. (2008). Planetary 
isochrones for various metal content, age and irradiation 
are also plotted on Fig. 14. As Fortney et al. (2007) 
note, planetary radii and their evolution are scarcely af- 
fected when irradiation is small, and a equiv > 0.1 AU, 
thus it is expected that HAT-P-llb will match non- 
irradiated models better. And indeed, the best match 
is with extreme metal rich non-irradiated Baraffe mod- 
els with Z — 0.9 fraction of heavy elements (rock and 
ice), and a 10% H/He envelope. The non-irradiated 
metal-rich models of various ages between 3 to 7 Gyr lie 
very close to each other, and are all within the obser- 
vational error-bars. However, the non-irradiated models 
with significant (50%) H/Hc envelope (top of figure) are 
far from matching the observations. Baraffe et al. (2008) 
notes that in general, the distribution of the heavy ele- 
ments (core vs. envelope) can have considerable impact 
on the planetary radius, but in this extreme metal rich 
scenario, the freedom to vary their distribution is limited. 
Altogether, HAT-P-llb appears to be a super-Neptune 
planet with Z — 0.9 and with a 10% H/He envelope. 

It is interesting to compare HAT-P-llb to GJ 436b 
- the only other super-Neptune with known mass and 



radius. These two planets are very similar in some phys- 
ical properties, in spite of the very different environ- 
ment. There are numerous system parameter determi- 
nations for GJ 436b available in the literature (see com- 
pilation and Table All of Southworth 2008). Start- 
ing from the discovery of the transits by Gillon et al. 
(2007) (22. 6M®, 3.95 i?©), planetary masses range from 
M p = 22.26 M® to 24.8 M®, and planetary radii range 
from 3.95 R B to 4.9 (Bean et al. 2008). Some of these 
earlier estimates, including that of Bean et al. (2008) 
are plotted on Fig. 14. The RV amplitude of GJ 436 
(K = 18.34 ± 0.52 ms" 1 , T08) is larger and better deter- 
mined than HAT-P-llb {K = 11.6 ± 1.2ms -1 ), and the 
transits are deeper (6.5mmag vs. 4.2mmag), and have 
even been observed from space by Spitzer (Gillon et al. 
2007; Deming et al. 2007) and HST (Bean et al. 2008). 
The reason for the larger K value and transit depth of 
GJ 436, however, are due to the host star being a small 
mass and radius M-dwarf (M* w 0.45 M Q , R* w 0.46 R Q , 
T08 and Southworth 2008). Uncertainties in both the 
planetary mass and radius remain considerable because 
of the uncertainty in the stellar parameter determina- 
tion for M dwarf stars (T08). Altogether, the masses 
and radii of GJ 436b and HAT-P-llb are very close (to 
within 2-<r). GJ 436b has a slightly smaller and bet- 
ter determined mass. HAT-P-llb has a comparable and 
somewhat better determined radius. 

There are other similarities between the two plan- 
ets. Both orbit on short period mildly eccentric orbits, 
GJ 436b: e = 0.16±0.02 (Gillon et al. 2007) and HAT-P- 
llb: e = 0.198±0.046. The fact that they have not circu- 
larized yet is interesting, as the circularization timescale 
appears shorter than the age of the host stars (e.g. Mat- 
sumura et al. 2008, and references therein), although the 
stellar ages are ill determined due to the unevolved M and 
K dwarf host stars. As Matsumura et al. (2008) noted 
on the origins of eccentric close-in planets, some of these 
planets may have larger than expected Q specific dissi- 
pation functions (depending on the planetary structure, 
and characteristics of the tides). The fact that neither of 
these hot Neptunes has circularized yet may tell us some- 
thing about the formation and planetary structure of 
Neptunes. We note that the putative HAT-P-llc planet 
(§ 6.2) causing the long-term drift in the data appears to 
have a large semi-major axis, because the RV drift does 
not show any significant non-linear coefficient over the 
1 year timescale of the observations (§ 2.4). Thus, it is 
unlikely that it is causing eccentricity pumping. Whether 
the Kozai mechanism (Fabrycky & Tremaine 2007) plays 
an important role will be subject to further investigations 
when the long-term drift in the RVs is better character- 
ized. Both GJ 436b and HAT-P-llb have similar sur- 
face gravities, \ogg p = 3.11 ±0.04 (cgs, T08) for GJ 436b 
vs. loggp = 3.05±0.06 for HAT-P-llb. They also appear 
to have similar mean densities, p p = 1.69 ±012 gcm~ 3 
(T08) for GJ 436b and 1.33±0.20 gem" 3 for HAT-P-llb. 

It is surprising that in spite of these similarities (mass, 
radius, density, surface gravity, eccentricity), the two 
planets are in grossly different environments. Conse- 
quently, the GJ 436b— HAT-P-llb pair must yield tight 
constraints on theories that should reproduce their prop- 
erties in spite of the different environments. GJ 436b 
orbits an M dwarf with half the mass and radius of 
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HAT-P-11, one tenth of the luminosity (L* = 0.026 Lq 
vs. 0.26 ± 0.02 Lq), and colder effective temperature of 
T eff * = 3350 ± 300 K (T08) versus 4780 ± 50 K for HAT- 
P-llb. The smaller semi-major axis of GJ 436b compen- 
sates to some extent for the lower luminosity host star, 
but still, the irradiation they receive is grossly differ- 
ent (GJ 436b: 3.2 • 10 7 erg s" 1 cm" 2 , HAT-P-llb: 1.34 • 
10 s erg s _1 cm~ 2 , i.e. 4 times more, as time-integrated 
over the eccentric orbit). These lead to different expected 
equilibrium temperatures: 650±60K for GJ 436b (T08) 
vs. 878 ± 15 K, assuming complete heat redistribution. 
Also, the different stellar effective temperatures lead to 
different spectral distributions of the infalling flux. It will 
be interesting to compare the atmospheric properties of 
the two planets to see the effect of different integrated 
flux and spectral flux distributions on these otherwise 
similar planets. 

While both GJ 436b and HAT-P-llb have similar mass 
and radius, and inferred structure of Z = 0.9 heavy el- 
ement content with 10% H/He envelope, it is notewor- 
thy that the metallicity of the host stars, and presum- 
ably the environment these planets formed in, is dif- 
ferent; [Fc/H] = -0.03 ± 0.20 for GJ 436 (T08) and 
[Fe/H] = +0.31 ± 0.05 for HAT-P-11. 

From the observational point of view, HAT-P-11 is a 
bright star (V <~ 9.58), more than a magnitude brighter 
than GJ 436. The transits of HAT-P-llb are at a small 
impact parameter (b = 0.347^q ^g) as compared to the 
near-grazing transit of GJ 436b (b = 0.848). The to- 
tal duration of the transit for HAT-P-llb is about twice 
(0.0957±0.0012 days) that of GJ 436b (0.042 d). We can 
thus expect very high quality follow-up observations of 
HAT-P-llb in the future. 

6.2. A Second Planet in the System? 

The radial velocity of HAT-P-11 shows, in addition 
to the 4.8878 day period induced by HAT-P-llb a sig- 
nificant trend of Gi = 0.0297 ± 0.0050 ms' 1 day- 1 . 
This drift could be due to a second planet in the sys- 
tem, corresponding to M 2 sin«2/ 'a\ ~ G\/G = (0.061 ± 
0.01)Mj kU~ 2 , where the subscript "2" refers to HAT- 
P-llc, and G is the gravitational constant. Such a con- 
jecture is scarcely surprising: for example Bouchy et al. 
(2008) have noted that 16 out of the 20 hot Neptune plan- 
ets discovered heretofore are members of multiple planet 
systems. Interestingly enough, as pointed out by Bouchy 
et al. (2008), this fraction (80%) is significantly higher 
than 23%, which is the fraction of all known exoplanets 
which are in multiple systems. We also recall that the 
significant orbital eccentricity of HAT-P-llb is reminis- 
cent of that of G J 436b, and that the latter has motivated 
the speculation that there may be another planet in the 
system, sufficient to maintain the eccentricity of GJ436b 
Ribas et al. (2008); Bean & Seifahrt (2008). Finally, as 
was also pointed out by Bouchy et al. (2008), of 7 multi- 
ple planet systems containing both a gas-giant planet and 
a hot Neptune planet, all have super-solar metallicities. 
The high metallicity of HAT-P-11, [Fe/H] = +0.31±0.05, 
then strengthens the suspicion that there may be multi- 
ple planets in the system. 

Continued monitoring of the radial velocity of HAT-P- 
11 over the next few years should indicate whether the 
drift persists, and whether it shows evidence of curvature 



of the orbit. Another way to establish the presence of a 
second body in the system is to search for transit timing 
variations in repeated transits (Agol et al. 2005; Holman 
& Murray 2005). 

6.3. Future Kepler Observations 

The fact that HAT-P-11 will lie on one of the detec- 
tor chips for the Kepler spacecraft presents remarkable 
opportunities for scientific followup of the discovery pre- 
sented here. If all goes well, Kepler should be capable of 
1-min cadence photometry at a precision of 0.1 mmag for 
this bright star. Currently Kepler is expected to be oper- 
ational for a minimum of 3.5 years, and perhaps longer; in 
that time it should have more than 250 transits by HAT- 
P-llb. In addition it should show photometric variability 
due to surface inhomogeneities such as starspots, span- 
ning more than 40 rotations if indeed the rotation period 
of the planet is verified to be about 29.2 days. Here we 
list some opportunities for important scientific followup 
with Kepler: 

Kepler transit data — With a depth of 4.2 mmag, the 
transit light curve of HAT-P-llb should be measurable 
by Kepler with excellent precision for a single light curve 
and extraordinary precision when many light curves are 
combined. To maximize the accuracy of parameters de- 
rived for the HAT-P-11 system, it is important that, 
in addition to the Kepler observations, precise ground- 
based radial velocity data be obtained through the du- 
ration of the Kepler observations. 

The possibility of a second planet in the system can 
be tested through a search for transit timing variations 
(TTV's) in the large number of transits to be observed 
by Kepler. As noted above, radial velocity data obtained 
over the same time frame could provide independent ev- 
idence for such a planet, and comparison of the two data 
sets would be of obvious importance. 

Kepler Out-of-transit photometry — The Kepler data 
should reveal the detailed modulation of photosphcric 
brightness-presumably due to starspot activity-with 
great precision and time coverage over the lifetime of the 
mission. This could yield much useful information about 
the star, including its rotation and even differential rota- 
tion (Frohlich 2007). Furthermore, if multiple transits of 
HAT-P-llb pass in front of the same star-spot on the face 
of the slowly rotating star, then extraordinary precision 
can be achieved in determining the rotation rate (Silva- 
Valio 2008). Comparison with concurrent ground-based 
measurements of the chromospheric emission S index will 
enable investigation of how photospheric spots correlate 
with chromospheric emission on stars other than the sun. 
As noted earlier, the data on S index reported here show 
a long-term (~ 450 day) variation reminiscent of solar 
activity cycles; monitoring both photospheric emission 
variations from Kepler broadband photometry and chro- 
mospheric emission variations from ground based tele- 
scopes could prove useful. 

Kepler asteroseismology — Kepler will obtain extraordi- 
nary information about the internal structure and evo- 
lution of stars from asteroseismic investigations, which 
measure the frequencies and amplitudes of acoustic os- 
cillation modes ("p- modes") through small changes in 
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Fig. 14. — The mass-radius diagram for Neptune-mass planets. Plotted are Uranus and Neptune from our own Solar system, the 
detection presented in this paper, HAT-P-llb, and the only other known transiting hot Neptune GJ 436b. For GJ 436b we plotted all 
available mass and radius determinations (Southworth 2008), the most recent one being the highest point in radius. Also shown in the 
upper right corner is the hot Saturn HD 149026b. Ovcrplottcd on the figure are theoretical models. HAT-P-llb is clearly much bigger 
than a ~ 25.8 ± 2.9 Mq rocky or water-rich "super-earth" would be, as based on models from Valencia, Sassclov, & O'Connell (2007). 
HAT-P-llb is also inconsistent with irradiated extreme metal-rich (Z = 0.9) models of Baraffe et al. (2008) (labeled as "Baraffe Z=0.9, 
I, 3,5,7 Gyr", "I" stands for irradiated). Indeed, these models assume a equ i v = 0.045 AU semi-major axis around a solar twin, instead of 
the low irradiation received by HAT-P-llb(a e9I1 i„ = 0.1 AU). The Z = 0.5 models from Baraffe ct al. (2008) with a 50% H/He envelope lie 
on the top of the figure, clearly presenting a mismatch with HAT-P-llb. However, HAT-P-llb is fully consistent with Z = 0.9 metal-rich 
non-irradiated ("Nl") Baraffe et al. (2008) models (thick lines), in particular with the 3-5 Gyr models. 

p-modes (e.g. Woodard & Noyes 1985; Woodard et al. 
1991), because the sub-surface structure is modified in re- 
gions of solar magnetic activity. Monitoring the p-mode 
frequencies of HAT-P-11 over the course of the Kepler 
mission, along with continued ground-based monitoring 
of the S-index, may show whether similar behavior occurs 
in the sub-surface layers of another star with significantly 
different structure. 

Finally, we note that two other stars with known tran- 
siting planets, namely TrES-2 (O'Donovan et al. 2006) 
and HAT-P-7 (Pal et al. 2008a), also lie on the Kepler de- 
tector array, so it should be possible to carry out studies 
of these stars and their planets similar to those described 
above for the HAT-P-11 system. It would be interesting 
to compare the results for these two stars, with spec- 
tral type F and G respectively, with those for the K star 
HAT-P-11. Almost certainly many more TEP systems 
will be discovered by Kepler, but those systems that are 
known in advance of the mission should be especially 
valuable since observations with the 1-minute cadence 
can be scheduled from the start of the mission, plus the 
longer baseline improves our sensitivity to detect secular 
changes. 



the integrated brightness of the star. The "large split- 
ting" of p-mode frequencies yields direct information on 
the radius of the star, and their "small splitting" gives 
information on the H/He ratio in the core and hence 
the age of the star (e.g. Christensen-Dalsgaard 2004). A 
great deal of additional information is possible if the star 
is accompanied by a transiting planet (like HAT-P-llb): 
in this case the transit light curves yield completely inde- 
pendent information on the radius and mean density of 
the star. An additional bonus in the case of HAT-P-11 is 
that it has an excellent Hipparcos parallax. Combining 
these independent pieces of information should lead to 
very detailed knowledge of its internal structure, includ- 
ing second-order near-surface effects not accounted for in 
normal interior structure models (Kjeldsen et al 2008). 
Information on the age of the star, if it can be obtained 
from measurement of the small splitting, could provide a 
useful comparison with the age inferred either from the 
global (light curve, RV, stellar isochrone) modeling of the 
data or from the Ca HK-index. 

As noted above, there is a suggestion from the varia- 
tion of its S index that HAT-P-11 is undergoing long- 
term variations of its surface magnetic activity some- 
what reminiscent of the solar activity cycle. On the Sun, 
solar cycle-related variations of solar magnetic activity 
are known to be correlated with frequency variations of 
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APPENDIX 



SYSTEMATIC VARIATIONS 



As regards the red-noise or systematic variations (also referred to as "trends"), we may distinguish between those 
effects that have a well understood reason, and where deviations correlate with a set of external parameters, and other 
systematic variations where the underlying parameters are not known. The deviations are usually taken with respect 
to the median of the light curve. 

The m(ti) light curve of a star observed at U time instances can be decomposed as: 



where mo(j>, U) is the systematics-free light curve as a function of time, and parametrized by a discrete set of parameters 
denoted as p, G{U) is white noise, E is the systematic variation due to the change in external parameters el, and T 
denotes the additional trend that is not known to be simple function of such parameters. For planetary transits the 
functional form of mo can be a transit model by Mandcl & Agol (2002), and the p parameters can be the a/R+ 
geometric semi-major axis, p = R+/ R p relative diameter of the planet, the normalized impact parameter b and the T c 
center of transit. For a 5 Scuti variable star, these parameters can be the amplitudes and phases of various Fourier 
components. 



Although briefly described earlier (Bakos et al. 2007b), it is worth defining the External Parameter Decorrelation 
(EPD) technique, since it is extensively used in this work. The EPD method attempts to determine the actual functional 
form of E in Eq. Al. The EPD effects are treated and determined as specific for each star, i.e. no information from 
other stars are used (this is a key difference when compared to TEA that uses a template of other stars — see later) . 

Constant (or simple) EPD: — The simplest form of EPD is when we assume that the underlying mo signal is constant, 
and also that (E) > (T) (i.e. the effects to be corrected by EPD are of the same order as those corrected by TFA). A 
typical application is for planetary transits in survey mode, since such transits are short events compared to the total 
orbital time, and have mostly been observed around stars with much smaller variation than the transit amplitude (note 
that there are exceptions; see Alonso et al. 2008), thus the underlying signal can be approximated with its median 

^95% of the time. The E relation between Sm(ti) = m(ti) — mo(i») = m(ij) — (m(tj)) m and ej(ij) parameters are 
sought via a chi-square minimization, usually by non-linear least squares method with outlier rejection ((),„ denotes 
median value). Upon determining E, we get an EPD corrected signal SmEPD(U) = Sm(ti) — E(ti). 

Our experience with HATNet survey data has shown that the effects of E and T are indeed comparable, and thus 
we use constant EPD. Specifically, the external parameters are the X,Y sub-pixel position, the background and its 
standard deviation, the stellar profile parameters characterizing the PSF width (called S) and its elongation (D and 
K), plus the hour-angle and zenith distance. 

Simultaneous EPD: — Used primarily in the analysis of the photometric follow-up data, where most of the observations 
are centered on the transits, thus the assumption of constant signal does not hold. One possibility is to use the 
out-of-transit (OOT) section of the light curve to determine the EPD parameters in constant EPD mode, and then 
apply the correction for the in-transit section of the light curve as well. This is sub-optimal, as the OOT may be very 
short compared to the in-transit section, or may be missing. Thus, the key difference compared to constant EPD is 
that we use all data to recover the dependence on external parameters, i.e. mo(p,t») in Eq. Al is not constant, but 
assumed to be a function of other parameters and time. The fitting procedure determines both the E correlation and 
the optimal set of p parameters simultaneously. 



The signal after (or without) the EPD procedure still contains the general systematic variations denoted as T{ti) in 
Eq. Al. These are suppressed by the Trend Filtering Algorithm (TFA, see Kovacs et al. 2005), assuming that certain 
other stars in an M element light curve "template" show similar variations. 

To recall, for selected star j, TFA minimizes the following expression: 



m{ti) = m {p, h) + G(U) + E{el{ti)) + T(U), 



(Al) 



External Parameter Decorrelation 



The Trend Filtering Algorithm 




(A2) 



where 



k=M 




(A3) 



k=l,k=£j 
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is the TFA filter function of M elements. The notation is like above in Eq. Al, mjfi(p,ti) is the model function and 
Ej is the EPD correction. 

There are many variants of TFA, and they have been used widely in this work. The signal search in the HATNet 
data is done in two parallel steps. First via simple (or constant) TFA where there is no assumption on the periodicity 
the model function rrijfi is a constant, Ej is zero, as EPD has been performed as an independent step before TFA, 
and the Ck coefficients are sought. Second, via reconstructive TFA, where the simple TFA is followed by a frequency 
search, the signal is phase-folded with the most significant frequency, the model function nij.o is fitted to the folded 
data, the model function is un-wrapped to the original time-base, and D is minimized again to iteratively determine 
the Ck coefficients. The third method, introduced in this work, is simultaneous TFA, whereby the TFA coefficients 
and the functional dependence of m,j(p,ti) on p parameters are simultaneously determined. 

EPD and TFA can be performed sequentially, or even simultaneously. In the analysis of the photometric follow-up 
data (see § 2.3) we implemented such a simultaneous EPD-TFA, where the Ej(e) EPD function is simultaneously 
fitted with the Ck TFA template coefficients and the rnj model function. 

Both TFA and EPD can be performed globally, using one set of coefficients for the entire dataset, and locally, when 
Eq. A2 is split up into smaller data-blocks, such as one-night segments, and the c& TFA coefficients and the EPD 
function parameters are fitted for each data-block separately to allow for changing systematics. Conversely for global 
TFA systematic variations of template stars have to "match" those of the main target (minus the model) for all nights. 

Finally we note that if we do not use EPD (eliminate the term for E in Eq. Al), but apply TFA only then naturally, 
TFA will take care of some of the systematics that otherwise would have been corrected by EPD. 



